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ABSTRACT 


Pulsed weather radars can be used to depict meteorological conditions such as 
rainfall rate and wind velocity. These quantities are calculated from measurements of 
reflectivity, mean radial velocity and velocity spread using echo signal samples from 
weather targets. These radar measurements derive from modified radio frequency (RF) 
echo signal parameters, including amplitude, frequency and phase, returned to the radar 
from the weather target. RF scattering and propagation effects modify echo signal 
parameters. Bias and variance in the weather signal parameter estimates naturally 
influence the accuracy of all subsequent quantities produced. For meteorological 
products to be as accurate as possible, the amount of uncertainty in each estimated 
quantity must be minimized. If radar system parameters are not accurately known, the 
reflectivity estimate will be biased. A well-controlled calibration process is therefore 
critical to reduce the bias of the reflectivity estimate. This thesis presents the design and 
implementation of one such calibration system, specifically for use with the MWR-05XP 
(a Mobile Phased-Array Pulse-Doppler X-band Weather Radar first created at the Naval 
Postgraduate School in 2005), although the general results are applicable to all radars. 
The calibration system presented is an active, external calibrator intended to verify end- 
to-end radar system performance. Within this thesis, a background of the radar 
calibration problem along with the research objectives for this specific project and related 
works are introduced and discussed. The theory behind the operation of weather radar 
(how the three principle quantities are measured and related to signal parameters) is also 
presented and explained. The density function for precipitation, relation between signal 
correlation and velocity spread, and fundamentals of weather radar signal parameter 
estimation are given. The specific design is described in explicit detail and measured 
results provided in appendices. Conclusions on the success of the implementation are 
drawn and recommendations for future work are presented. 
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EXECUTIVE SUMMARY 


Pulsed weather radars can be used to quantify meteorology into amounts such as 
rainfall rate maps and wind velocity fields. These quantities are calculated from 
measurements of reflectivity, mean radial velocity and velocity spread using echo signal 
samples from weather targets. These radar measurements derive from modified radio 
frequency (RF) echo signal parameters, including amplitude, frequency and phase, 
returned to the radar from the radar’s target. RF scattering and propagation effects are 
the mechanisms which modify echo signal parameters. Bias and variance in the weather 
signal parameter estimates naturally influence the accuracy of all subsequent quantities 
produced. For meteorological products to be as accurate as possible, the amount of 
uncertainty in each estimated quantity must be minimized. If radar system parameters are 
not accurately known, the reflectivity estimate will be biased. A well-controlled 
calibration process is therefore critical to reduce the bias of the reflectivity estimate. This 
thesis presents the design and implementation of one such calibration system, specifically 
for use with the MWR-05XP (a Mobile Phased-Array Pulse-Doppler X-band Weather 
Radar first created at the Naval Postgraduate School in 2005), although the general 
results are applicable to all radars. The calibration system presented is an active, external 
calibrator intended to verify end-to-end radar system performance. 

The approach is a repeater that will produce a known Radar Cross Section (RCS) 
via loop gain control and create Doppler velocities and Doppler spread via phase 
modulation with a sawtooth. Proper control of the sawtooth frequency will also allow 
introduction of Doppler spread. Although a constant amplitude return would suffice, use 
of a programmable, constant phase attenuator would also allow control of the amplitude 
statistics of the return to the radar. Thus, with full control of both amplitude and phase, 
one could replicate the return for each of the Swerling target types, although it is only 
precipitation return that is really of interest in our application. A block diagram of the 
implementation of this proposed Active Calibrator System (ACS) design is presented 
below. 
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The reflectivity for a weather radar contains radar system parameters which may 
not be precisely known. The objective of calibration is to determine the composite of 
these parameters. The composite of these factors is also known as the radar system 
calibration constant. In Dr. Jeffrey Knorr’s paper, “Weather Radar Equation Correction 
for Frequency Agile and Phased Array Radars,” it has been shown that, for the case 
where the reference target is on boresight, 6 = 0, the measurement is made at the 
reference frequency, X = Xq, and the range is short {La = 1) the radar system calibration 
constant can be described as: 





4 J 


1 

b 

_1 


From this equation, it is clear that a measurement of receiver output power and 
range for a reference target of known cross-section {6 = 0, X = Xo, La = 1) permits the 
radar system constant on the left side to be determined and the radar calibrated for 
reflectivity estimation. Phase modulation is required to produce Doppler shift, according 
to: 


v{t) = Acos'\joj + O{t) =y4Re 




For 0(t) periodic where the phase modulated signal has a spectrum given by a sawtooth 
waveform of amplitude adjusted to produce a phase change that is a multiple of In (see 
waveform pictured below) then, 

v(t) = where (A6’ = m2;r). 
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The signal is translated (Doppler shifted) from frequency to frequency . 



By cycling through 360° (In) of phase at this periodic rate a known 
frequency/Doppler shift is produced against which the weather radar can be calibrated. 
This periodic phase cycling (and thus frequency translation) is implemented in the ACS 
using a digital phase shifter. Velocity spread can be added to the Doppler shifted signal 
by varying the frequency of the sawtooth waveform used in the phase modulator. The 
frequency should be varied in accordance with a Gaussian random variable that has been 
low pass filtered to produce a signal with the desired correlation time. 

On the basis of these principles, the ACS was constructed (RF box internal 
configuration digital photograph follows): 



All control of the ACS was achieved through National Instruments Lab VIEW 
(user interface shown below). 
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After construction, measurements were made to determine the equivalent 
reflectivity (RCS) of the ACS. The result of those measurements and that entire testing 


process are presented below. 
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Time vs. Amplitude for Empty Chamber, Sphere and ACS. (ESC_plotTA) 
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Empty Chamber, Sphere, & Calibrator - Frequency Domain 



GHz 
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Frequency vs. Amplitude for Empty Chamber, Sphere and ACS. (ESC_plotF) 

The first plot is a composite of all data sets in the time domain, not zoomed (all). 
The second plot is a composite of all data sets in the frequency domain. All data sets 
include the measured data for the three targets: an Empty Anechoic Chamber (chamber), 
the 12” Aluminum Sphere (sphere) and the ACS (calibrator). Timing gates were placed 
around the sphere’s location (42 ns) and calibrator’s location (57 ns). These timing gates 
eliminated undesirable echoes. One such undesirable echo was at <2 ns. This echo 
corresponds to the coupling between transmit and receive antennas which were placed 
about 18” apart. There are seven graphed data sets on each plot: the sphere without a 
time gate and with a time gate centered on its 42 ns return, the calibrator without a time 
gate and with a time gate centered on its 57 ns return, the empty chamber without a time 
gate and with a time gate for the calibrator, then with a time gate for the sphere. 


Eor the sphere, the time gate (G2) was centered at 41.8 ns and spanned 39.8 to 
43.8 ns. Eor the calibrator, the time gate (Gl) was centered at 56.8 ns and spanned 54.8 
to 58.8 ns. The time gates correlate directly to target return distances. The distance 
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(range) from transmit and receive antennas in the anechoic chamber to the target pedestal 
inside the chamber is 19 feet (or 38.6 ns in time). The distance is traversed twice (once to 
the target and then from the target) at the speed of light. The reason the time gate for the 
sphere was observed at and consequently set at 42 ns is the sphere was suspended from 
the ceiling of the chamber above and slightly behind the target pedestal inside the 
anechoic chamber. This added approximately an additional 1.6 feet of distance and 3 ns 
of time because the hypotenuse rather than the longest leg of the triangle was traversed. 
The reason the time gate for the calibrator was observed at and consequently set at 57 ns 
is different than for the sphere. This time gate was approximately 8 ns greater because of 
the time delay due to transmission line lengths rather than target distance. The 8 ns is 
almost 8 feet (one-way path), slightly greater than the straight line distance between 
transmit and receive antennas of the ACS (as expected). See the example calculations 
below for the mathematical relation between time and distance (two-way path). 

c = 3.00x10* m/s 

i? = 5.7912 m (19 feet) 

9 /? 

— = 38.6 ns 
c 

t = 41.8 ns 

^ = 6.2700 m (20.6 feet) 

Initially, calculations were made to predict the received power from both the ACS 
and 12” Sphere. The expected (calculated) value for received power from the ACS was 
-36.5dBm. The measured value for received power from the ACS was -34.0dBm. The 
expected (calculated) value for received power from the 12” Sphere was -47.9dBm. The 
measured value for received power from the 12” Sphere was -46.7dBm. The expected 
(calculated) values very nearly equaled the measured values. 

The next plot is from the data in the frequency domain of both the sphere and 
calibrator previously presented, now used to calculate the equivalent RCS of the 
calibrator. Measurements were made from 9-10 GHz. The Vector Network Analyzer 
used in the test setup as the signal source and receiver in the anechoic chamber recorded 
the magnitude of the returned amplitude power in logarithmic form (log MAG) in units of 
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milliwatts (dBm). All measurements are in the frequeney domain gated for time (FG). 
This ineludes both the sphere (SFG) and ealibrator (CFG). The measured values were 
very close to the expected (calculated) values. The received power for the 12” 
Aluminum Sphere was calculated to be -47.8679 dBm and observed to be near -47 dBm. 
The received power for the calibrator was calculated to be -36.5032 dBm and observed to 
be near -35 dBm. The RCS of the sphere was previously calculated to be -11.3688 dBsm 
(S-RCS). To calculate the RCS of the calibrator (C-RCS) based on the 
observed/measured data, the CFG - SFG=A (dB) was added to the sphere’s RCS. Over 
the entire frequency band the C-RCS varied from 0.6 to 3.2 dBsm. Note all these 
measurements in the chamber included the use of the RF amplifier within the calibrator 
system. When used with the weather radar, signal levels are great enough and ranges 
short enough such that an RF amplifier is not needed. Furthermore, calculations are 
made to make sure the signal return to the radar is in the proper range to prevent receiver 
saturation. 


Equivalent RCS of Calibrator 
(with in-line amplifier) 



Frequency (GHz) 


With these measured values the equivalent RCS of the ACS was observed to 
range from 0.6 to 3.2 dBsm across the 9-10 GFlz RF band. On average (mean), the 
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equivalent RCS of the ACS was 1.4934 dBsm. The median was 1.5101 dBsm. The high 
value was 3.1742 dBsm at 9 GHz and the low value was 0.5862 dBsm at 9.38 GHz. 

An aetive ealibrator system (ACS) for weather radar was sueeessfully designed 
and implemented. The ACS was verified to be a funetional weather radar aetive 
ealibration system under laboratory testing. It sueeessfully provided the means (met the 
requirements) to verify all three of the signal parameters measured (estimated) by a 
weather radar. These inelude the zeroth Doppler moment (refieetivity), first Doppler 
moment (average radial veloeity), and seeond Doppler moment (veloeity spread). 
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I. 


INTRODUCTION 


A. BACKGROUND 

Pulsed weather radars produee estimates of refleetivity, mean radial velocity and 
velocity spread using echo signal samples from a wide variety of weather targets. From 
these radar return signals, other meteorological quantities such as rainfall rate and wind 
field characteristics are also derived. Estimates are derived from the parameters of the 
modified radio frequency (RF) echo signal (amplitude and frequency/phase) scattered 
from the target back in the direction of the radar. RF scattering and propagation effects 
are the mechanisms which modify radar echo signal parameters and it is those 
modifications that contain the desired weather-related information as a modulation of the 
carrier wave. Bias and variance in the weather signal parameter estimates naturally 
influence the accuracy of all subsequent quantities produced, and condition assessments 
made. For meteorological products to reflect weather conditions as accurately as 
possible, the uncertainty of each estimate must be minimized, or accurately characterized. 
If radar system parameter accuracies are not precisely known, the reflectivity estimate 
will be biased and thus, assessments of true conditions could be either misleading or 
masked. A well-controlled calibration process is therefore critical to reduce reflectivity 
bias due to imprecise knowledge of the radar system parameters. There may also be bias 
associated with the statistical estimation process and its resulting estimator, but poor 
calibration is the most serious source of error. The variance of the estimate of interest 
will depend on the number of samples and their correlation. 

B. RESEARCH OBJECTIVE 

The research objective of this thesis is to design and implement a radar calibration 
system (calibrator). The calibrator is for specific use with the Naval Postgraduate 
School’s MWR-05XPi X-band weather radar, but the general concept is applicable to all 

1 World’s first and currently (2007) only Mobile Weather Radar created in 2005 operating at X-band 
in Pulse-Doppler with a Phased-Array. In acronym, this is the MWR-05XP . 
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radars. The calibrator presented is an active, external calibrator intended to verify end-to- 
end radar system performance. An ideal active calibration system produces a radar return 
with known RCS and Doppler shift permitting the determination of the radar calibration 
constant. As of 2005, the MWR-05XP operates with a calibration constant determined 
through the use of radar system parameters that have previously been only approximately 
known. A known RCS return condition can be obtained when a metal sphere is used to 
calibrate a radar. This calibration method involves the suspension and isolation of the 
sphere from great height at long distance^ in ideal weather conditions. These pristine 
conditions are not always available when and where the MWR-05XP is in use. In fact, 
the MWR-05XP is primarily located at the Center for Interdisciplinary Remotely Piloted 
Aircraft Studies (CIRPAS) in Marina, CA on the grounds of the municipal airport (a 
relatively small workspace) where the heights and distances required for calibration are 
never available. The calibration of the MWR-05XP is done periodically on a large test 
range at Camp Roberts, CA over 100 miles away. The calibrator research involved in 
this thesis effort will eliminate the prohibitive burdens (physical conditions and 
workload) required with flying the sphere while simultaneously allowing more specific 
calibration parameter determinations along with the inherent advantage of more frequent 
calibration constant verification. The calibrator also allows the added calibration of 
Doppler-related quantities (which is not currently possible). The final advantage is the 
end-to-end method. It uses the entire radar system in normal operating mode. 
Calibration systems sometimes rely on piecemeal measurements that are combined 
together for performance estimates, thus leading to larger uncertainties. This system does 
not. 

Research in this effort will implement an existing design of an active calibrator 
for the MWR-05XP and will not involve the creation of an innovative design 
configuration. It will neither involve design optimization, nor investigation of 
alternatives. This calibrator, once implemented, will be used with the MWR-05XP 
exclusively (i.e. no other radars involved) to determine the radar calibration constant. 

2 Sphere heights of 0-1000m with ground distance equal to 700m have been used with the MWR-05XP 
to cover beam elevation angles from 0-55°. 
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The type of signal used by the calibrator in this process will be of the minimum required 
complexity, although the performance parameters involved are those that most accurately 
reflect operational weather characteristics. 

The primary hardware involved with this thesis effort will be a digital phase 
shifter, programmable attenuator, and real-time controller with reconfigurable chassis. 
The only software involved will be Lab VIEW (in the manipulation of the controller 
hardware) and MATLAB (in the analysis of data and display of results). 

C. RELATED WORK 

The thesis by Schroder presents a weather radar signal simulator. According to 
the abstract of the Schroder thesis, 

the simulator was developed in MATLAB and implemented several 
different functionalities allowing for stepped frequency, multiple pulse 
repetition frequencies (PRFs), pulse compression using a chirp waveform, 
and variation of both weather and radar input parameters. Post-detection 
processing capabilities include autocorrelation and Fast Fourier Transform 
(FFT - for single PRF only); estimation of weather parameters such as 
reflectivity (Z); average Doppler, radial velocity, and velocity spread; 
pedagogical plots including a phasor plot of phase change over time and a 
velocity histogram, instantaneous observed reflectivity and power for each 
pulse over time.[l] 

The output of Schroder’s simulator was modified to provide the calibrator controller 
software with Doppler frequencies and amplitude levels indicative of real-world weather 
effects for use in the calibration process. 

The article by Marshall and Hitschfeld investigates the random fluctuations of the 
radar echo produced from a random array of scatterers (such as weather scatterers). 
According to the introduction to this article, “the intent is to establish (theoretically and 
by analogue computation) the extent of the limitation these fluctuations have on the 
accuracy of the measurement information obtainable. Additionally, what improvements 
may be effected by averaging or by modifying the method of pulsing or scanning are 
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explored.”[2] The Marshall and Hitschfeld theory is used in this researeh to produee the 
amplitude levels indieative of real-world weather effeets for use in the ealibration 
proeess. 

The report by Pratte and lee explores external radar ealibration methods. 
Aeeording to the exeeutive summary of the Pratte and lee report, “the targeted radar for 
this report is the WSR-88D, the weather radar in use (as of 2007) by the U.S. National 
Weather Serviee in its national weather radar network. It reviews availability and 
feasibility of “ealibration aids”, equipment and software that will assist in verifying 
WSR-88D refleetivity measurement and in redueing its uneertainty.”[3] The 
reeommendations of the Pratte and lee report were used to enhanee understanding of the 
general weather radar ealibration problem and its alternative solutions. 

The text by Doviak and Zrine emphasizes the applieation of Doppler radar for the 
observations of stormy and elear air. Aeeording to the introduetion (1.2) to the text, “it 
presents a eomprehensive treatment of the teehniques used in extraeting meteorologieal 
information from weather eehoes and relates radar and signal eharaeteristies to 
meteorologieal patterns.”[4] The Doviak and Zrine text served as a elassie referenee 
book on radar theory and aeeepted teehniques applied to meteorology. 

The paper by Knorr presents a derivation of two eorreetion faetors to the elassie 
Probert-Jones weather radar equation for use with frequeney agile, phased array weather 
radars. Aeeording to the abstract of this paper, “the result is an extended weather radar 
equation for use with these advanced weather radars. The two derived factors account for 
the effect of a changing RF carrier and electronic beaming steering not present in classic 
weather radars when computing a reflectivity estimate from weather signal data 
samples.”[5] The Knorr paper, as well as personal discussions with its author (an advisor 
to this thesis), provided the necessary understanding of the unique calibration 
requirements of these advanced weather radars from theory to practical application. 

Chapter Four of the Federal Meteorological Handbook discusses potential 
operational uses of the WSR-88D products as well as the structure of certain 
meteorological phenomena. According to the introduction to Chapter 4 (4.1), “it 
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provides information on observing and foreeasting meteorologieal phenomena using this 
radar’s products.”[6] Of particular interest in the development of this thesis’s research 
was the discussion in early sections of operationally troublesome radar returns due to 
such things as ground clutter and sidelobes. The pictorial representations in the 
Meteorological Handbook were invaluable. 

D. ORGANIZATION OF THE THESIS 

There are four main chapters and several appendices which comprise this thesis. 
Chapter I gives a background of the radar calibration problem along with the research 
objectives for this specific project. Related works are introduced and discussed. An 
executive summary is also included. Chapter II explains the theory behind the operation 
of weather radar (how the three principle quantities of reflectivity, velocity and velocity 
spread are measured and related to signal parameters). The density function for 
precipitation, relation between signal correlation and velocity spread, and fundamentals 
of weather radar signal parameter estimation are given. Chapter III describes the specific 
calibrator design in explicit detail. Measured results are presented and explained in 
appendices. Chapter IV draws conclusions on the success of the implementation and 
provides recommendations for future work. 
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II. BACKGROUND AND SUPPORTING THEORY 


This chapter presents and explains the theory behind the operation of weather 
radar. The diseussion includes how the three prineiple quantities of reflectivity3, mean 
radial veloeity^ and veloeity spread^ are measured and are related to radar signal 
parameters. The density function for precipitation, relation between signal eorrelation 
and veloeity spread, and fundamentals of weather radar signal parameter estimation are 
given. 

Significant portions have been exeerpted from Chapter 2 “Weather Radar and 
Weather/Rain Theory” of Ulf Schroder’s thesis “Development of a Weather Radar Signal 
Simulator to Examine Sampling Rates and Scanning Schemes.” Both Schroder’s thesis 
and this thesis are in support of the larger researeh project/interest of Dr. Jeffrey Knorr 
(advisor to both theses). Much of the background and supporting theory used in 
Schroder’s thesis (eompleted 2005) applies to this thesis (completed 2008) and the larger 
researeh project/thrust. 

A. WEATHER RADAR VS. TRADITIONAL RADAR 

Weather radar samples the scattered RF energy from distributed targets such as 
precipitation, hail, snow, melting iee, and clear air. This is unlike traditional radar which 
views these returns as clutter, unwanted eehoes that interfere with the observation of the 
return from a point target of interest. Traditional radar usually implements Doppler 
filtering to diseriminate between slow moving weather clutter and fast moving point 
targets. In this application, the scattered energy contains descriptive characteristics of the 
weather condition. 


3 Function of the scatterer intensity. 

Function of the mean component of scatterer motion in the radial direction from the radar. 

3 Function of the variability of Doppler velocity values due to turbulence within and velocity shear 
across the sampling volume. 
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B. RADAR RESOLUTION CELL 


Weather radar differs from traditional radar in another important way. Traditional 
radar seeks to encapsulate its target of interest (approximately a point target) entirely 
within a single radar resolution cell. For weather radar this is neither practical nor 
informative at either the micro or macro scales. It is impractical to image a single rain 
drop or an entire storm and even if it were possible the information provided would not 
be especially useful. Instead, weather radars image a volumetric segment of the 
meteorological phenomena of interest. This volume is defined by a segment of 
propagated energy (see Figure 1) covering pulse width, x. The depth of this segment 
depends on the pulse width according to: 



Figure 1. Illustration of Radar Resolution Cell. [From 1] 

Schroder’s figure above illustrates and later expressions (specifically, reflectivity 
factor Z) use the volume of the resolution cell. It is uniquely defined here for 
convenience, as: 
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( 2 . 2 ) 


AK = n=£(^«,)(^A)(f)45) 

See Chapter II, Seetion E, Subseetion I, “Antenna Gain and Beam effeetive Solid 

Angle,” for an explanation of the last faetor, ——. 

2 In (2) 

C. SWERLING TARGET TYPES 

Meteorologieal eonditions produee a radar target return that varies with time. 
Besides the target itself, fluetuations in the target’s return may also be eaused by 
variations in the medium between the target and the radar (atmospherie and 
meteorologieal eonditions not the intended target), radar system instabilities (platform 
motion and equipment instabilities), target aspeet ehanges, and many other variables. In 
order to analyze all of this variation, the target is modeled as an overall system. For 
systems analysis purposes, only the ‘gross’ behavior of a target needs to be known, not 
every detail of all the complicated physics involved in the complex scattering situation. 
Therefore, fluctuations caused by all of these many sources (exclusive of those caused 
simply by a change in range) are usually grouped together and represented by a single 
quantity for a specific target type. 

The following excerpt was relayed from the radar text by Skolnik, although nearly 
every radar text contains some explanation of Swerling target types. Let cr be a random 
variable with a probability density function (PDF) that depends on the factors associated 
with target fluctuations other than a change in range and the overbar represents the 
average value of the random variable. Two PDFs are commonly used; 

p{cr) = Le-^l^ (2.3) 

cr 

and 

p{cT) = , PDF #2 (2.4) 

cr 

These two PDFs approximate the nature of the scattering of RF energy from many 
targets. There are targets which consist of many independent scattering elements of 
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which no single one (or few) dominates (PDF #1) and targets whieh have one main 
seattering element that dominates in the presenee of smaller independent seattering 
sourees (PDF #2). These two PDFs along with the rate at whieh the fluetuations take 
plaee allow the definition of four types of radar targets, also known as Swerling target 
types. Swerling target types I and II are assoeiated with PDF #1 and types III and IV 
with PDF #2. The rate of the fluetuations in this seattering is the other differentiating 
faetor whieh further defines the Swerling target type. There are targets with slow 
fluetuations (sean-to-sean) and targets with rapid fluetuations (pulse-to-pulse). Swerling 
target types I and III are of the former and types II and IV of the latter. 


Swerling Target Types and Probability Density Functions 


PDF#1 

PDF #2 

Slow Fluetuations 
(sean-to-sean) 

Type I 

Type III 

Rapid Fluetuations 
(pulse-to-pulse) 

Type II 

Type IV 


Table 1. Swerling Target Types and Probability Density Funetions. 

As a weather radar, the MWR-05XP targets Swerling types I and II. The 
ealibrator deseribed in this thesis was eonsequently developed for those same Swerling 
types, but the general proeess of the design and implementation of a radar ealibrator ean 
be applied to any radar and Swerling target not just weather radars and Swerling types I 
and II. 

D. SAMPLE CORRELATION 

To paraphrase Sehroder, the power level of seattered RF energy from distributed 
targets varies over time. An estimate of the average power level is aehieved by averaging 
power over a number samples. The varianee of a refleetivity estimate ean be improved 
simply by averaging more samples. Flowever, these samples may not be independent of 
one another. Independenee signifies a suffieient level of ehange in the distribution of the 
target within the volume being imaged. For this ‘suffieient’ level of ehange to have 
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occurred there exists a minimum required period of time. This minimum time period 
guarantees samples are not eorrelated from one observation period to the next. 


Schroder explains, “The power spectral density of a meteorological signal is 
approximately Gaussian and ean be written as 


jf-ff 

2al 


(2.5) 


S{f) = S,e 

where f is the Doppler frequency and cr^ is the Doppler spectrum variance. Taking the 
Fourier transform of Equation (2.5) and normalizing results in the correlation coefficient 


p{t)i 




where is the time spectrum variance and is related to cr^ by 

1 


cr, = 


2n(7 1 


and 


^/ = 


2cr„ 




( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 


By using the correlation coefficient a measure of minimum difference in time can 
be estimated ensuring dependent or independent samples. Nathanson elaims that for 
independent sampling /7(r)<0.02 is required whereas /7(r)>0.15 will ensure 

dependenee. Doviak and Zrnic set a correlation threshold for coherence at /7^(r)>e ' 
which corresponds to a value of /7(r) > 0.6 . Using this threshold renders 

Z 


T; <0.25- 


2rcr„ 


(2.9) 


For highly eorrelated samples." 

A level of ‘sufficient’ ehange has oeeurred within an observation period for a 
sample so that any subsequent samples may be deemed independent if: 

Z 


T >0.70- 


nG„ 


( 2 . 10 ) 
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where is the sample time period and <j^ is the RMS veloeity spread. 

E. ESTIMATION OF POWER DENSITY AND REFLECTIVITY (Z) 


I. Antenna Gain and Beam effective Solid Angle 


Antenna gain varies depending on azimuth and elevation angles for direetive 
antennas. For a point target, antenna gain is eonstant only for a fixed azimuth and 
elevation angle defined between the radar and a point in spaee at a speeific time. For a 
volumetrie target, antenna gain varies aeross a range of azimuth and elevation angles 
defined between the radar and a volume of spaee at a speeifie time. The antenna pattern 
funetion, f{ 6 ,<j)), integrated over the solid angle eorresponding to the radar resolution 
eell is used to estimate a value of antenna gain for volumetrie targets. 


From Sehroder (and Probert-Jones as the original souree), “The magnitude of the 
antenna pattern funetion is first squared to yield a power relation and squared again to 
agree with the two way propagation of radar. 

\\f{e,<l>)\dCL ( 2 . 11 ) 

Att 


The Gaussian Pattern Model is the antenna pattern funetion that will be used. It 
provides a good approximation of the main beam between the half power points. 
Extending it in both azimuth and elevation provides 


= \e\<U, \<i>\<U 


( 2 . 12 ) 


where = —— and 


41n2 


41n2 


and 6*j, are the half power beam widths in the 


orthogonal planes. 


Performing integration over the solid angle gives 

+00 +00 

\\f{ 04 )\dn=\\e 


le^ 

1 




Ak 


ded(t> = ^^{\-A) 

81n2^ ’ 


(2.13) 
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where ^ < 0.034 and represents the eontributions from the sidelobes. In a further 
approximation, at its maximum value the (l-^) term eontributes a negligible 0.15 dB 
and therefore ean be ignored.” 


The final form of the expression used to approximate beam solid angle for 
volumetrie targets is 




4;t 


dQ. = 


n6^(j)^ 

81n2 


(2.14) 


See also the Probert-Jones paper[7]. The equation above gives beam solid angle 
eontributing to the volume return. Note this result is based on a Gaussian approximation 
to the beam pattern. As Dr. Knorr shows in his paper, integration of an aetual antenna 
pattern will generally result in a correetion whieh may not be negligible. 


2. Physics of Precipitation 

The following presents some basies in preeipitation physies relevant to 
developing a physieal understanding of preeipitation, speoifieally rain the primary target 
of weather radar. The following three seetions have been seleetively ehosen from a 
larger diseussion of the physies of preeipitation by Sehroder. Eaeh of the three seetions’ 
text has been paraphrased and equations reprodueed verbatim. 


a. Size and Shape 

Raindrop size determines the terminal veloeity at whieh the drop falls and 
the shape of the drop. Small drops, with a diameter D < 0.35 mm, are spherieal. Medium 
size drops, with diameters 0.35 mm < D < 4 mm, have progressively flattened bases and 
are approximated as spheroids. Large drops, D > 4 mm, do not last very long and tend to 
break apart upon eollision. The largest drops, D > 10 mm, are unstable and break apart 
without eollision. 
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b. Drop Size Distribution (DSD) 


Radar cross section (RCS) density, r/, is the average RCS density per unit 
volume, as previously mentioned, and for independent seattering is the sum of the RCS of 
all seatters in the return volume divided by that volume. The DSD represents the 
diameter density (mm/m ) of drops for all diameters. The Marshall-Palmer DSD 
provides a general model: 

N{D) = N,e-^ 

A = 4.1R “"W ' (2.15) 

= 8 X10^ m^^mm ' 

Applying truneation with a maximum drop diameter of yields 


/V(D) 


\N,e-^,D<D^^ 

lo, 


(2.16) 


c. Rainfall Rate 


Rainfall rate determines the depth of aeeumulated water per unit time. It 

ean be derived from water eontent and fall speeds. A eloud’s water density is 

00 

M = (2.17) 

® 0 

where is the water density of a rain drop. 

Turning this into rainfall rate 

00 

R = -\D^N[D)v[D)dD. (2.18) 

Rainfall rate is estimated from refleetivity using any one of a number of 
slightly different models, Z = Ar^. Z = 300r' for example, is the eurrent WSR-88D 
default model. The best rainfall model depends upon the type of storm or equivalently, 
drop size distribution. Dual polarization radars provide more information about drop size 
and thus rainfall rate. Latest results from JPOL experiment (dual polarization WSR-88D 
at NSSL) are given in Ryzhkov[8]. 
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3. 


Weather Radar Range Equation and Estimation of Reflectivity 


Again, following Schroder, the basie form of the radar range equation (negleeting 
minor losses): 


R = 


{47rf R" 


(2.19) 


“In the basie form of the weather radar equation radar eross seetion (RCS), a, is 
replaeed by RCS density, rj , whieh is the expeeted RCS density per unit volume. Sinee 
the target is a volume target, the expeeted return needs to be integrated over the whole 
volume, whieh ean be divided into the solid angle and depth of the return. Letting the 
power level at the target be 

and integrating the return power over the return volume renders 


p,=\dp,. 


(4^ 


ffujl/lejpa 

R ^ Q 


1024;z-" R^\n2 


( 2 . 21 ) 


reeognizing that R ? 


CT 

~ 2 ' 


Radar meteorologists use the reflectivity factor, Z, instead of RCS density where 

7 = ^K,rz (2.22) 

and 

(2-23) 

^ ieV 

Substituting Equation (2.22) into Equation (2.21) and solving for Z yields the refleetivity 
estimator 


1024;z-^RMn2 Z" 


— 1024RMn2Z" 

^ P,G,G^CT0,M'Kf ' 


(2.24) 
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For an assumed continuous distribution of drops sizes with diameter density 
N{D) Equation (2.23) becomes 

^max 

Z = I D^N{D)dD. (2.25) 

0 

Solving the integral form for an assumed Marshall-Palmer Drop Size Distribution 
(DSD) Equation (2.15), yields 

00 ^ I JLT 

Z={D^N(D)dD = -(2.26) 

i ^ ^ (4.ii?--')' 

Drops do not exist in sizes ranging up to infinite diameter, but are likely to be 
limited by a maximum diameter. Doviak and Zrnic state that Equation (2.26) 
overestimates Z and suggest a truncated DSD with D = . In this case. 


Z = 




(4.1i? 


-r{7,a) 


(2.27) 


where a = AD^^^^and y{l,a) is the incomplete Gamma function. Equations (2.26) and 


(2.27) provide an analytical means for computing Z.” 


F. RADAR PRINCIPLES AND WEATHER RADAR EQUATION 
CORRECTION 


The basic results for phased array weather radar are presented from Knorr’s 
paper. This paper addresses issues unique to frequency hopping, phased array weather 
radars and quantifies the need for a target of known radar cross section (RCS). 

“The first term in the equation below indicates those radar system parameters 
which may not be precisely known. 

~ PP^Gl 

4 

where 

Xo = reference^ wavelength (m) 


cos^ 


cr 


= R. 


(2.28) 


^ The choice of reference wavelength within the operating band of the radar is arbitrary. 
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Ls = radar system losses (ratio) 

La = atmospherie propagation loss (ratio) 

Pout ^ output power (W) 

R = range to target (m) 

Pt = transmit power 
Grx = ree elver gain 

Go = gain at the referenee frequency,/o 
L = operating wavelength (m) 

6 = angle of the antenna beam with respect to boresight or the normal to the array face 
a = target radar cross section (RCS) (m^) 

The objective of calibration is to determine the composite of these factors with a 
single measurement. To simplify this process, several conditions can be selected. These 
include: calibration on boresight with the antennas of the active calibrator system (ACS) 
and the radar (0 = 0), calibration at the reference RF center frequency (^ = ^o), and 
calibration at a range that is short enough that atmospheric losses can be ignored (La = 1). 
Then, the equation simplifies and rearranges to: 

~ P,G^Gl 

4 

With this relationship, a measure of the radar receiver’s output power and range 
for a reference target of known cross-section (the ACS’s equivalent reflectivity) will 
determine the radar system constant (left hand side of the equation).” 


(4;r) 


' R^P.. 




(2.29) 


For frequency agile, phased array weather radars, the equation for reflectivity 
takes the form 


10241n2Zo 

1 

4 

L P 

out 

7r^{cT)\K^A^ 


_PtG^Gl 

_(///o)"cos^_ 


where 

Z = reflectivity (mm /m ) 

Ao = reference wavelength (m) 

Ls = radar system losses (ratio) 

La = atmospheric propagation loss (ratio) 

P = average output power (W) 

(JUl 

R = range to target (m) 
c = speed of light in a vacuum = 3x10^ m/s 
X = pulsewidth (s) 


(2.30) 
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r 


r 


.93 (H^O) 

Sr = hydrometeor relative permittivity 

^„„and (3„„are the broadside 3 dB prineipal plane beamwidths at the referenee 

B\J £>\j 

frequeney 

Pt = transmit power 

Grx = ree elver gain 

Go = gain at the referenee frequeney,/o 
/= operating frequeney (Hz) 
fo = referenee’^ frequeney (Hz) 

6 = angle of the antenna beam with respect to boresight or the normal to the array face 

Note that the units of Z are given correctly (mm /m ) by the first two terms of 
Equation (2.30) if range is in millimeters, Zg in the numerator of the first term is in 
millimeters and the quantities in the denominator of the first term are in meters. 

This is the corrected form of the classical Probert-Jones classical weather radar 
equation as derived by Dr. Jeffrey Knorr in his paper, “Weather Radar Equation 
Correction for Erequency Agile and Phased Array Radars.” 

The third term contains radar system parameters which may not be precisely 
known. The objective of calibration is to determine the composite effect of these factors. 
The composite of these factors is also known as the radar system calibration constant. In 
Dr. Knorr’s paper it has been shown that, for the case where the reference target is on 
boresight, 6 = 0, the measurement is made at the reference frequency, 2 = Xq, and the 
range is short {La= 1) 

4 

Prom Equation (2.31), it is clear that a measurement of receiver output power and 
range for a reference target of known cross-section {6 = 0, X = Xq, La = 1) permits the 
radar system constant on the left side to be determined and the radar calibrated for 
reflectivity estimation. 


( 4 ^)'^' 4 . 

Xlcj 


(2.31) 


^ The choice of reference frequency within the operating band of the radar is arbitrary. 



G. 


ESTIMATION OF VELOCITY 


According to Schroder, “Pulsed radar ean, by measuring the phase ehange 
between eonseeutive pulses, estimate radial veloeity. The effeet, ealled Doppler, is the 
frequeney ehange that oeeurs due to the relative movement of a target with respeet to the 
radar. Considering a two-way propagation path renders a total phase ehange of 

2 R 

(f> = 2n — (2.32) 

X 


In terms of angular frequeney the result is 

_d(l) _ An dR _ Anv^ 
dt X dt X 




(2.33) 


where is the radial veloeity and is the Doppler frequeney. 


From a pulse pair, using a eoherent deteetor, the phase ehange ean be extraeted 


using 

V = I + jO = A 

n n J n 

^^Z{Vn]=<l>n 

V V*=A ^ 

^ n+V n 

leading to the single pulse pair veloeity estimate 


4;r 


(2.34) 


(2.35) 


where T is time between pulses and / is pulse repetition frequeney.’ 


Also from Sehroder, eoherenee between pulse pairs is required to measure phase 
ehange over time for a volumetrie target to minimize phase ehange from effeets other 
than radial movement. The varianee of a veloeity estimate ean be improved simply by 
averaging more samples. Building on the single pulse pair veloeity estimate Equation 
(2.35), the multiple pulse pair estimator beeomes: 


$ 

V = - 


1 N-\ 

— 

iv-is 


1 


N-\ 


N-\ An tt 




(2.36) 
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H. ESTIMATION OF VELOCITY SPREAD 


Turbulence, wind shear, and differences in fall speed due to drop size are but a 
few of the many directly influencing factors on the spread in radial velocity of 
precipitation. There are even indirect factors (motion of the radar antenna) which impart 
an apparent spread in velocity. The net effect of these factors is the width of the spread 
in the velocity spectrum and the time for de-correlation between samples. 


Schroder explains, “An estimate of the velocity spread can be computed using the 
velocity samples obtained from a pulse pair comparison between single pulse pair 
velocity estimates Equation (2.35). The estimator for sample variance is given by 




m 


i=\ 




N 


and the estimator for standard deviation is” 

. 

1. CALIBRATION REQUIREMENTS 


(2.37) 


(2.38) 


Factors exist both internally (hardware and software) and externally 
(environmental and remote sensing effects) that influence a radar’s measurements apart 
from the target return that is trying to be quantified. The calibration process is used to 
reduce the influence of these factors on the composite measurement. Internal calibrators 
minimize system variations by separating the instrument from variable external 
influences (such as through the use of known reference signals generated internally to test 
internalized system components). External calibrators, however, operate completely 
outside (and independent of) the instrument under test. By operating completely outside, 
external calibrators incorporate those elements in the calibration process that describes 
external influences that would effect the instrument’s measurements, but are not 
considered with an internal calibration method. By operating independently, external 
calibrators are able to calibrate for the same factors targeted by internal calibrators 
without the calibrator being influenced in its calibration by the instrument under test. 
With these two additional traits, external calibrators can verify whole-system calibration. 
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The result is a true “end-to-end” ealibration of the specifie radar eonfiguration for the 
speeifle external measurements that are about to be aequired. 

Beeause these internally and externally influeneing faetors are so many and so 
varied they eannot be individually quantified in theoretieally expressions. Instead, these 
eombined effeet faetors are often denoted by, in terms of Pratte, “the so-ealled “radar 
eonstant” C. In any given implementation, the radar eonstant C ean eontain slowly 
varying quantities, such as receiver conversion gain, system and site losses (to include 
propagation path loss caused by the atmosphere), and “hidden” assumptions. In WSR- 
88D software jargon the radar constant derived from the engineering calibration is called 
“SYSCAL” (Sirmans 1992), and temporal variations obtained by the internal, electronic 
calibration-maintenance function software is called “delta SYSCAL” or “CAL#”.” 

J. CALIBRATION METHODS AND SELECTED METHOD 

A summary of alternative calibration methods was obtained from Pratte. A list of 
these methods follows: 

1. Traceability 

2. Solar and Celestial Radiative Flux Transfer Methods 

3. Radar Calibration Targets 

4. Lunar Reflectivity 

5. Transponder Methods 

6. Elevated Ground-range Methods 

7. Near-Field Measurements 

8. Reflectivity Comparison Methods 

9. Rain gauge Comparison Methods 

10. Dual Polarization Consistency Methods 

Compared to other calibration methods, use of an external active calibration 
system is the most reliable, independent and thorough means of calibration. For this 
reason, this type of system was chosen for implementation. 

The approach is a repeater that will produce a known RCS via loop gain control 
and add Doppler return information via phase modulation with a controlled sawtooth. 
Proper control of the sawtooth frequency introduces Doppler spread, an important 
measurement quantity in the weather radar. This approach assumes a smooth profile in 
Doppler return as a function of the varying sawtooth. A calibration sphere is selected to 
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enable that smooth response profile (see Figure 6). Although a eonstant amplitude return 
would seemingly suffiee, use of a programmable, constant phase attenuator would also 
allow software control of the amplitude statistics of the return to the radar. Thus, with 
full control of both amplitude and phase, one could replicate the return for each of the 
Swerling target types, although it is only the precipitation return that is really of primary 
interest in our application. 

K. SUMMARY 

This chapter has lead up to the formulations of the radar equation for a point 
target and for a volume target which show that the radar calibration constant can be 
determined if, having a target of known RCS, the echo power from that target is 
measured. The objective of this thesis’s research is to determine if an active calibrator 
can be built which will produce the desired known RCS over a response range of interest 
for weather returns. Introducing a known Doppler shift and known RMS frequency 
spread allows one to compare the radar estimates of those parameters with the values set 
for the calibrator. The Doppler shift also addresses problems associated with use of the 
calibrator in an environment where there is fixed clutter return. 
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III. DESIGN, IMPLEMENTATION AND OPERATION 


A. DESIGN 

This study implemented an existing design (see Figure 2) of an active calibrator 
system (ACS) for the MWR-05XP. It did not involve the creation of such a design nor 
did it involve the investigation of alternatives to or optimization of this design. 



Figure 2. Block Diagram of the ACS. 

See Figure 3 for a layout of the actual components used. These components 
comprise the internals of the Active Calibrator System or ACS that was implemented in 
this study. 
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The Active Calibrator (internal) 
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Figure 3. Diagram of Component Layout inside the ACS Enclosure. 


B. HARDWARE IMPLEMENTATION 

1. System Description Overview 

See ACS Pictures Appendix. 

The active calibrator system (ACS) is comprised of two Pasternack sectoral horn 
antennas mounted on either end of a 6 ft long l”x4” piece of lumber. To this piece of 
lumber, another piece of lumber of the same dimensions is attached to create an “L” 
bracket. Midway between the two antennas is centered the ACS enclosure, an aluminum 
box (17”xl3”x5”), to house the power supplies, control computer, and RF components 
(programmable attenuator, digital phase shifter, etc.). On the face of this enclosure, 
parallel to the antennas’ mainbeams, is attached a second piece of lumber (cut to the 
length of the enclosure face) to present a flush, all wood box exterior. Along this exterior 
and running the entire length of the “L” bracket between the two antennas, X-band Radar 
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Absorbent Material (RAM) is attached. Although this RAM has its greatest RF 
absorbency at 9.4 GHz, it provides sufficient RF absorbency across the entire RF range of 
interest, 9-10 GHz, to reduce undesired RF returns greater in power than the transmitted 
power of the ACS’s transmit antenna output. This RAM is sectioned into three pieces 
(two 24”xl0”x2” and one 19”x9”x2”). Each piece is attached via removable Velcro to 
aid in transport and storage. The antennas are mounted and oriented to transmit and 
receive horizontally polarized signals. They are connected to the enclosure via 3 ft RF 
cables. All of the ACS’s components apart from the antennas are contained inside the 
enclosure. 

2. Hardware List 

External Hardware (Quantity/Description/Dimensions, if applicable) 

(1) Aluminum box(enclosure) (17”xl3”x5”) 

(3) Pieces Emerson & Cuming AN77 RAM (two, left and right side, 24”xl0”x2” and 
one, center, 19”x9”x2”) 

(3) Pieces lumber (l”x4”, two 6 ft long and one 19” long) 

(2) Pasternack Sectoral Horn Antennas 

(2) Pasternack semi-flexible RE Cables (3 ft) 

(1) Ethernet (RJ-45) Cable 

(1) Dell Precision M65 with 2.33 GHz Core 2 Duo Processor, 2GB 667 MHz RAM, 80 
GB HDD Laptop and Dell 90W-AC Power Supply 
(1) AC Power Cable 

Miscellaneous Mounting Hardware (nuts, screws, washers, bolts, etc.) 

Internal Hardware (Quantity/Description) 

(1) Laptop Power Supply (24V CUI EA1050B-240) 

(1) Condor DC Power Supply (GLC-50A) 

(2) 25-pin-to-15-pin Control/Power Cables 
(1) Voltage Regulator Circuit 

(1) National Instruments Control Computer (CompactRIO Model Line to include; (1) 
real-time controller, (1) reconfigurable embedded chassis, and (2) I/O modules) 

3. Voltage Regulator Circuit 

A voltage regulator circuit was required to convert the DC power supplied by the 
Condor power supply (-15VDC, +15VDC, and +5VDC) to the additional voltages (- 
5VDC and -I-12VDC) required by the RE components. 
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Voltage Regulator Circuit Hardware List (Quantity/Description) 
(2) 0.01 uF Capacitors 
(1)33 ohm Resistor 

(1) 7905C (provides -5VDC from -15VDC input) 

(1) 7812CT (provides +12VDC from +15VDC input) 



Figure 4. Voltage Regulator Circuit Diagram. 

4. RF Components 

See Appendix on Digital Phase Shifter. 

See Appendix on Programmable Attenuator. 

See Appendix on Antennas. Appendix describes antenna specifications. 

See Appendix on Antenna Testing. Appendix describes antenna patterns, gain 
and beamwidth calculations, and isolation. 

5. Control Hardware 

Control of the ACS is achieved through the use of the National Instruments (NI) 
Compact Reconfigurable Input/Output (CompactRIO) product line. From the 
ni.com/compactrio website, “The National Instruments CompactRIO programmable 
automation controller (PAC) is a low-cost reconfigurable control and acquisition system 
designed for applications that require high performance and reliability. The system 
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combines an open embedded architecture with small size, extreme ruggedness, and hot- 
swappable industrial I/O modules. NI CompaetRIO is powered by reconfigurable I/O 
(RIO) FPGA technology.” Again from the ni.eom/pac website, “PACs eombine 
programmable logie controller (PLC) ruggedness with PC functionality under an open, 
flexible software arehitecture.” 

The CompaetRIO products used by the ACS include: a real-time controller, a 
reconfigurable embedded ehassis, and two I/O modules for digital input. Software 
development was done on a Windows PC laptop with NTs Developer Suite with Real¬ 
time deployment and FPGA deployment options installed. Through this developer suite 
and these options, all software development was done in the NI Lab VIEW™ visual 
programming language. For items not run on either the laptop’s Windows 0/S or the 
real-time controller’s real-time 0/S, LabVIEW was auto-complied into the VHSIC 
Hardware Description Eanguage (VHDE) and run on the reconfigurable embedded 
chassis. This process was seamless to the developer and no specific knowledge of VHDE 
was required. These items included all of the software control required to manipulate the 
phase shifter and programmable attenuator on a time scale at or near its fastest possible 
switehing speed. A normal 0/S or real-time 0/S is incapable of this speed of control 
requiring the use of Eield Programmable Gate Arrays (EPGAs). 

See Appendix on CompaetRIO. 

6. Hardware Limitations 

There are a few hardware limitations associated with the ACS. These limitations 
affect the maximum possible Doppler shift and therefore radial velocity and velocity 
spread that can be practieally generated. Eimited factors are dietated by the switehing 
speed of the programmable attenuator (300 ns) and phase shifter (500 ns), the CRIO 
EPGA clock (25ns), and the response time of the digital I/O module (100 ns). Obviously, 
the switching speeds of the programmable attenuator and phase shifter are the true 
limiting factors as the CRIO components are more than eapable of switching states at or 
faster than the fastest switching speeds of both RE components. 
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See the example calculation of the effect of the phase shifter’s fastest possible 
switching speed. 

1 


2'x(500xl0 


= 7812.5 Hz Doppler Frequency Shift 


299792458 

7812.5 X-= 117.2 m/s Radial Velocity 

2x9990000000 

117.2 X 2.23693629 = 262.2 mph Radial Velocity 
This limitation actually falls below what is physically possible (or what 
has been observed) and, therefore, could potentially limit ACS applicability in a few 
applications. The fastest wind speed ever recorded on Earth was 318 mph in one of the 
tornadoes that hit the suburbs of Oklahoma City on 3 May 1999.^ However, the fastest 
recorded wind speed not associated with a tornado was measured at 230 mph in 1934.9 


ACS FREQUENCY/DOPPLER CALIBRATION AND SOFTWARE 
IMPLEMENTATION 


In order to calibrate the weather radar’s ability to measure radial velocity and 
velocity spread, the ACS must be able to produce a known Doppler shift and Doppler 
spread. By comparing the known quantity to the measured quantity and adjusting for any 
discrepancy, the calibration can be completed and the weather radar’s accuracy assured. 
This section describes how the ACS achieves this. 


The active calibrator’s loop gain determines equivalent reflectivity, Zeq.^*’ Phase 
modulation is required to produce Doppler shift, according to: 


v(t) = Hcos[(»^t + 6’(t) =HRe 




(3.1) 


Following the derivations in Tan’s thesis [9] and the paper by Cummings [10], for 
0(t) periodic. 


^ National Weather Serviee Foreeast Offiee, "Frequently Asked Questions about the May 3, 1999 
Bridge Creek/OKC Area Tornado," retrieved 2008-09-21 from 
http://www.srh.noaa.gOv/oun/storms/19990503/may3faqs.php 

9 University Corporation for Atmospherie Researeh (UCAR), "Windows to the Universe/Wind," 
retrieved 2008-09-21 from http://www.windows.uear.edu/tour/linkAearth/Atmosphere/wind.html 

See seetion on “ACS Amplitude/Refleetivity Calibration” (next seetion) for a diseussion on the 
determination of equivalent refleetivity. 
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Figure 5. Two Cyeles of the Periodie, Linear, Modulating Sawtooth Function. 
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sin(A6’/2-n;r) 


(3.4) 


[A0 / I-htt) 

In general, the phase modulated signal has a general signal given by 


(A6’/2-n;r) 


(3.5) 


If the amplitude of the sawtooth is adjusted to produce a phase change that is a multiple 
of 271 then 

A0 = mlTT for (m e /) 

rifor77 = m . (3.6) 

T [O for 77 + m 

Therefore, 

v(z) = where (A^ = 7?72;z-). (3.7) 
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The real part of the eomplex exponential expression is cosine with argument indicating 
the translation (Doppler shift) from frequency to frequency + mf^ . In other words, 
each spectral line is separated by the modulating frequency and lines exist symmetrically 
above and below the carrier. 


By cycling through 360° (2n) of phase at this periodic rate a known 
frequency/Doppler shift is produced against which the weather radar can be calibrated. 
This periodic phase cycling (and thus frequency translation) is implemented in the ACS 
using a digital phase shifter. However, there are imperfections associated with this 
method. Imperfect implementation results in spurious sidebands due to finite sawtooth 
flyback time, incorrect sawtooth amplitude, the digital staircase approximation to linear 
phase shift, phase error, and finite phase shifter switching time. [10] 


Again, following the derivations in the thesis by Tan [8], for a digital phase shifter 


using N=2^ steps. 


/ \ (“!)”* sin(;r/A) . . 


(3.8) 


^^mN + \ [n/N) 

where K = mN +1 and m is an integer. The frequencies of the non-zero spectral 
components are given by / = /^ +-^„, or f = +[mN + \) . The amplitude of the K* 


sideband relative to the input carrier is given by 


P^_ \C,fA^/2 , 

P, A^/2 ' 


The first sideband in the spectrum of the translated signal is 




(') 


sin(;r/ a) 
(;r/ a) 


Hcos(®^ + . 


The translation loss is defined as 


(3.9) 


(3.10) 


TL,s=-20log 


sin(;z-/ A) 
(;r/ a) 


= -20 log 


sin(^;r/2^ j 
(;r/2") 


(3.11) 


This loss rapidly approaches 0 dB as the number of bits of phase shift increase. If B=4 
bits, for example, TLdB=0.06 dB. 
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The amplitudes of the undesired sidebands are also of interest. The undesired 
sideband with the highest amplitude is K = -N + 1. Relative to the desired sideband, 
K=l, the suppression ratio is 

=201og(iV-l) = 201og(2''-l). (3.12) 

The suppression ratio also improves rapidly as the number of bits is inereased. If N= 4, 
for example, SRdB= 23.5 dB. 

See Appendix I on Lab VIEW for the Digital Phase Shifter whieh gives further 
explanation of how the phase shifter is eontrolled by the CompactRIO to produce the In 
periodic phase cycling. 

In conclusion, one of the primary advantages of the ACS is its ability to introduce 
Doppler shift. This permits the calibrator signal to be extracted from mainbeam ground 
clutter at elevation angles near zero and from sidelobe clutter at higher elevation angles. 
In most cases, convenience will dictate that the ACS be operated near the ground. 

D. ADDING VELOCITY SPREAD 

Velocity spread can be added to the Doppler shifted signal by varying the 
frequency of the sawtooth waveform used in the phase modulator. The frequency should 
be varied in accordance with a Gaussian random variable that has been low pass filtered 
to produce a modulating signal with the desired correlation time. This functionality was 
not fully implemented in the ACS at the conclusion of this study. See Chapter IV on 
Future Work. 

E. ACS AMPLITUDE/REFLECTIVITY CALIBRATION 

In order to calibrate the weather radar’s ability to measure reflectivity, the ACS 
must be able to reliably simulate a known reflectivity. This is accomplished by receiving 
the incident radar signal and retransmitting a signal of known power relative to the 
incident signal. From a measurement of the echo signal power received by the radar and 
knowledge of the ACS loop gain, the calibration can be completed and the weather 
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radar’s absolute accuracy in amplitude assured. This section describes how the ACS is 
calibrated so it can in turn be used to calibrate a weather radar. 

1. Measurement Premise 


The ACS needs to produce a power density equivalent to that scattered by a target 
of known RCS. This power will be the effective radiated power (or ERP) of the ACS. 
The ERP of the ACS will be a function of (depend on) the ERP of the radar. This is 
because the ACS does not produce an ERP from an independent/internal power source, 
but rather as acts a repeater (and modifier) of the power received from the radar. 


The power received from the radar will depend on the ERP of the radar and the 
losses involved between the ACS and the radar. At this radar frequency, the losses 
between the ACS and the radar will primarily depend on the distance between the ACS 
and the radar (spreading loss or free space loss only, unless a ground reflection is present 
as will be the case for elevation angles near zero) and atmospheric attenuation 
(propagation loss), but also losses internal to the radar (system losses which cannot be 
precisely known). Those losses which depend on distance will affect the reflectivity 
observed by the radar due to the calibrator twice - first, as the ERP of the radar is 
affected on its way to becoming the received power of the calibrator and, second, as the 
ERP of the calibrator is affected on its return to becoming the received power of the 


radar. 


While the spreading loss is well known {L = 



2 

) and can be accurately 


calculated, the propagation loss is less well known (as it depends on many time variant 
factors) and is more difficult to accurately calculate. In order to isolate this uncertainty at 
any given time, the equivalent reflectivity of the ACS can be compared to a known 
reflectivity at equal distances and with the propagation losses minimized as much as 
possible (as in an anechoic chamber). By comparing these two values, the equivalent 
reflectivity of the ACS can be measured. Under these settings the losses internal to the 
radar can be isolated (the system losses which cannot be precisely known will be 
calibrated out). 
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An equation for the ERP of the ACS and thereby the equivalent reflectivity of the 
ACS can be defined as a function of the received power of the ACS (from the radar) and 
the loop gain (system gain and losses) of the ACS. These include the gain of both 
antennas (transmit and receive), cable losses, phase shifter losses, programmable 
attenuator losses, isolator losses, amplifier gain, and other RF components’ losses (only 
the amplifier provides gain). On the whole, the ACS’s ERP and therefore its equivalent 
reflectivity can be controlled (varied) by the programmable attenuator (variable 
gain/loss). However, this is only effective as long as the net gain/loss of all of the ACS’s 
RE components (fixed gain/loss) in total exceed the reflectivity of the ACS due to its 
exposed surface area. 

The equivalent reflectivity due to the ACS’s transmit antenna is in direct 
competition with the actual reflectivity due to the ACS’s surface area. Both depend on 
the radar ERP, but the equivalent reflectivity further depends on the ACS’s loop gain 
acting to re-radiate the radar’s power while the actual reflectivity depends on the exposed 
surface area of the ACS (to include all sides facing the radar) acting to reflect the radar’s 
power. If the controlled power does not exceed this value, variable amplitude calibration 
(as will be needed for amplitude statistics) will not be possible. Instead, only a single 
reflectivity (the actual reflectivity) of the ACS will ever be observed by the radar. This is 
the same principle used in passive calibration systems (such as the balloon and sphere) 
where a controlled Doppler calibration is not possible. Radar Absorbent Material (RAM) 
was used to reduce the actual reflectivity to levels below the equivalent reflectivity. See 
Chapter III, Section B, Subsection I for further details on RAM use. 

2. Measurement Method 

What follows is the method used to calibrate out these unknown system losses for 
accurate reflectivity measurements (estimates) by the weather radar based on the 
equivalent RCS of the ACS. See Chapter II, Section F, on “Radar Principles and 
Weather Radar Equation Correction’’ for background and supporting theory directly 
applicable to this process and these calculations. The following development shows how 
the equivalent RCS of the ACS was experimentally determined through anechoic 
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chamber measurements and subsequent ealeulations to satisfy a most fundamental radar 
calibration requirement - namely, that the ACS’s equivalent RCS must be known [5]. 

As mentioned before, the value of the equivalent refleetivity (or more aeeurately 
equivalent radar eross seetion or RCS) of the ACS ean be measured and caleulated in 
eomparison to a known refleetivity (or RCS) of another target. In this study, an 
aluminum sphere was used and measurements made inside an aneehoie ehamber at a 
fixed, known distanee to serve as the known refleetivity. A 12” diameter aluminum 
sphere was ehosen out of three diameters available based on the initial ealeulations and 
estimates performed below. 


a. Initial Calculations 

Prior to making measurements in the aneehoie chamber theoretieal 
ealeulations were made using approximately known or individually measured values of 
the required quantities. The signal souree in the aneehoie chamber was a vector network 
analyzer (VNA). The VNA replaces the radar as the signal souree when ealibrating the 
ACS in the aneehoie chamber. 

Below is a list of eonstants, either measured or known, to include; 

R, the measured one-way distanee from the transmit and reeeive antennas to the ACS 
inside the aneehoie ehamber. (measured) 

f, the radio eenter frequeney. (known, seleeted on veetor network analyzer or VNA) 
c, the speed of light, (known) 

2, the wavelength, (caleulated) 

Pt, the VNA transmit power a.k.a. the power input to the aneehoie transmit antenna from 
the VNA to include all losses from the output of the VNA to the antenna (cabling, etc.), 
(known, selected on VNA) 

Gt=Gr, the gain of the VNA antenna - one in the same - a.k.a. the aneehoie transmit and 
reeeive antennas - both the same, (measured. See Appendiees on Antennas.) 

EIRP, the radar’s Effeetive Isotropic Radiated Power. The produet of the transmit power 
and antenna gain, (ealeulated) 

Gc, the ealibrator’s antenna gain. Both antennas on the ealibrator are the same, 
(measured. See Appendices E and E on Antennas.) 

Gamp, the ealibrator’s loop gain to include all system losses (eabling, phase shifter while 
not phase shifting, eonneetors, ete.) - i.e. net gain. The aetual gain of the amplifier used 
in the ealibrator is +35 dB with -20 dB of system losses for +15 dB net. (ealeulated via 
measured values; amplifier gain and system losses) 

FSL, the Tree Space Eoss or Spreading Eoss associated with the distance R. (calculated) 
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Cl, the calibrator’s input signal power as measured into the system at the receive 
antenna’s output port, (calculated) 

RI, the VNA input signal power as measured into the system at the receive antenna’s 
output port, (calculated) 

/ = 9.50xl0® Hz 
c = 3.00x10* m/s 

;L = —= 3.16x10^" m 

/ 

=25 dBm =316.2278 mW 
G, =G^ =16 dBm =39.8107 ratio 
EIRP = P,G, = 41 dBm = 12589.2541 mW 
= 13 dBm = 19.9526 ratio 
^amp =15 dBm = 31.6228 ratio 
7? = 5.7912 m (19 feet) 


2 

= 1.8830x10 ’ ratio =-67.2516 dB 


Cl = EIRP X FSL xG^= A.129% x 10^" mW = -13.2516 dBm 

RI = Cl X G„„,^ xG^x FSL x G, = 2.2371 x 10 " mW = -36.2032 dBm 


Constants 

Quantity: 

Range 

Frequency 

Speed of Light 

Wavelength 

Radar Transmit 

Radar Antenna 

Radar EIRP 

Calibrator 

Calibrator 






Power 

Gain 


Antenna Gain 

Amplifier Gain 

Symbol: 

R 

f 

c 

A=c/f 

p, 

G,=G, 

EIRP=P,*Gt 

Gc 

^amp 

Units: 

ft 

Hz 

m/s 

m 

dBm 

dB 

dBm 

dB 

dB 

Value: 

19 

9.50E+09 

3.00E+08 

3.16E-02 

25 

16 

41 

13 

15 

Alt. Units: 

m 




mW 

ratio 

mW 

ratio 

ratio 

Alt. Value: 

5.79 




316.23 

39.81 

12589.25 

19.95 

31.62 


Table 2. Summary of Constants Used in Estimated Equivalent RCS Calculations. 


FSL = 


A 


AnR 


Calculations 

Quantity: 

Symbol/Equation: 

Units: 

Vaiue: 

Free Space Loss 
FSL=[A/(4 *tt*R)]^ 
ratio dB 

1.88E-07 -67.25 

Caiibrator Input 
CI=EIRP*FSL*Gc 
mW dBm 

4.73E-02 -13.25 

Radar Input 
RI=CI*G3,^p*Ge*FSL*G, 
mW dBm 

2.24E-04 -36.50 


Table 3. Summary of Estimated Equivalent RCS Calculations. 


-36dBm is approximately the value of VNA input signal power expected 
when measuring the ACS inside the anechoic chamber. This will be the total returned 
power from the ACS due to the relayed power and not the reflected power. The reflected 
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power is much smaller (no return was detected with the ACS in the anechoic chamber 
and RAM attached) and has been minimized by the use of X-band Radar Absorbent 
Material (RAM). Again, see Chapter III, Section B, Subsection 1 for further details on 
RAM use. 

The power received from the ACS must be compared to the received 
power from a target of a known RCS to calculate an equivalent RCS for the ACS. The 
known RCS that was used was that of a polished aluminum sphere. Three spheres were 
available, a 3”, a 6” and a 12” diameter sphere. The expected received power was 
calculated for all three. The relation for the monostatic RCS of a sphere in the optical 
region (frequency region) was used. The power had to be great enough to be detected by 
the VNA, but not so great as to cause saturation. The best suitable candidate from the 
three available was found to be the 12” diameter sphere because of its larger (but not too 
large, i.e. saturating) return amplitude. See example calculation below, 
a = 6 in = 0.1524 m 

Pa = —^ = 30.3 ^ Optical Region (See Figure 6.) 

A. 

= 1 for Optical Region (See Figure 6.) 
na 

G = na^ =7.2966x10^ sm =-11.3688 dBsm 
Received Power = 1.6338x10^ mW = -47.8679 dBm 

Following the example, for a 3” diameter sphere -59.909 dBm and for a 6” diameter 
sphere -53.889 dBm. 
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Sphere Calculations 

Quantity: 

Sphere 

Radius 

Radius 

OPTICAL 

OPTICAL 

RCS 

RCS 

Received 

Received 


(Diameter) 







Power 

Power 

Symbol: 

- 

a 

a 

|3a=2TTa/A 

a/rra^ 

a = Tia^ 

a = Tia^ 



Units: 

- 

in 

m 

- 

- 

m^ 

dBsm 

mW 

dBm 

Value: 

3" 

1.5 

3.81 E-02 

7.58 

1 

4.56E-03 

-23.41 

1.02E-06 

-59.91 

Value: 

6" 

3 

7.62E-02 

15.16 

1 

1.82 E-02 

-17.39 

4.08E-06 

-53.89 

Value: 

12" 

6 

1.52E-01 

30.32 

1 

7.30 E-02 

-11.37 

1.63E-05 

-47.87 


Table 4. Summary for Calculations of Radar Cross Section of a Sphere. 


b. Anechoic Chamber Measurements 

Measurements were made in an anechoic chamber with the Vector 
Network Analyzer system. Three different series of measurements were made. One with 
the chamber empty (E), one with the sphere in the chamber (S), and one with the 
calibrator in the chamber (C). Measurements were made and represented both in the time 
domain (T) and frequency domain (F).ii To eliminate any possible interference from 
reflections (echoes) in the chamber due to less than ideal performance from the anechoic 
chamber (there were design compromises made in the chamber at the time of 


11 The basic measurement is a swept frequency (frequency domain) measurement. The time domain 
response is then obtained via TFT. Time domain resolution is determined by frequency domain bandwidth. 
A 1 GFlz BW should give about 1 ft. resolution. See reference [10]. It provides further explanation of this 
information. 
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construction due to space limitations) time based gating was used to eliminate undesired 
signals. Measurements were made both with this time based gating (G) and without (N). 
In the time domain, data sets and plots were often zoomed to the time/distanee of interest 
(when/where the calibrator’s return was expected). These are annotated zoomed (Z) or 
not (A, for all). When not zoomed, they inelude all data starting from time zero. Table 5 
summarizes the abbreviations, just discussed, that are used when referring to the various 
anechoic chamber measurement data sets. In all cases, a total 90 ns time base was used 
(RCS90). 


S=sphere 

T=time 

domain 

G=gated 

A=all 

(0 see to zoomed) 

Gl=54.8-58.8 ns 

E=empty 

chamber 

E=frequency 

domain 

N=not gated 

Z=zoomed 
(around time gate) 

G2=39.8-43.8 ns 

C=calibrator 



NOTE: does not apply 
to E, always 9-10 GHz 



Table 5. Naming Convention Summary. 


Figure 7 is a composite of all data sets in the time domain, not zoomed 
(all). Figure 8 is a eomposite of all data sets in the frequency domain. All data sets 
inelude the measured data for the three targets: an Empty Aneehoie Chamber (chamber), 
the 12” Aluminum Sphere (sphere) and the ACS (calibrator). Timing gates were placed 
around the sphere’s location (42 ns) and calibrator’s location (57 ns). These timing gates 
eliminated undesirable echoes. One such undesirable echo was at <2 ns. This echo 
corresponds to the coupling between transmit and receive antennas which were placed 
about 18” apart. There are seven graphed data sets on each plot: (1 and 2) the sphere 
without a time gate and with a time gate eentered on its 42 ns return, (3 and 4) the 
calibrator without a time gate and with a time gate centered on its 57 ns return, (5 and 6) 
the empty chamber without a time gate and with a time gate for the calibrator, then (7) 
with a time gate for the sphere. 
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Figure 7. Time vs. Amplitude for Empty Chamber, Sphere and ACS. 

(ESC_plotTA) 
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Empty Chamber, Sphere, & Calibrator - Frequency Domain 



GHz 


. Not Gated - S 

- Gated - S 

.Not Gated - C 

- Gated - C 

. Not Gated - E 

Gated for Calibrator - E 
- Gated for Sphere - E 


Figure 8. Frequency vs. Amplitude for Empty Chamber, Sphere and ACS. 

(ESC_plotF) 

Eigure 9 and Eigure 10 are data sets in the time domain of the sphere and 
the calibrator only, respectively. Each plot is of the respective target only and zoomed to 
very near the time gate size of each target. Eor the sphere, the time gate (G2) was 
centered at 41.8 ns and spanned 39.8 to 43.8 ns. Eor the calibrator, the time gate (Gl) 
was centered at 56.8 ns and spanned 54.8 to 58.8 ns. The time gates correlate directly to 
target return distances. As mentioned before, the distance (range) from the transmit and 
receive antennas in the anechoic chamber to the target pedestal inside the chamber is 19 
feet. 19 feet is 38.6 ns in time. The distance is traversed twice (once to the target and 
then from the target) at the speed of light. The reason the time gate for the sphere was 
observed at and consequently set at 42 ns is the sphere was suspended from the ceiling of 
the chamber above and slightly behind the target pedestal inside the anechoic chamber. 
This added approximately an additional 1.6 feet of distance and 3 ns of time because the 
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range hypotenuse of the distance (height geometry is involved) rather than the longest leg 
alone. The reason the time gate for the calibrator was observed at and consequently set at 
57 ns is different than for the sphere. This time gate was approximately 8 ns greater 
because of the time delay due to transmission line lengths rather than target distance. The 
8 ns is equivalent to 8 feet (one-way path), slightly greater than the straight line distance 
between the transmit and receive antennas of the ACS (as expected). See the example 
calculations below for the mathematical relation between time and distance (two-way 
path). 

c = 3.00x10® m/s 

i? = 5.7912 m (19 feet) 

9 /? 

— = 38.6 ns 
c 

t = 41.8 ns 

^ = 6.2700 m (20.6 feet) 
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Sphere - Time Domain 
(Aii=Zoomed) 



Not Gated 
Gated 


Figure 9. Time vs. Amplitude for Sphere, Zoomed. (S_plotTA=Z) 
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Calibrator - Time Domain 
(AII=Zoomed) 



Figure 10. Time vs. Amplitude for ACS, Zoomed. (C_plotTA=Z) 

The expected (calculated) value for received power from the ACS was - 
36.5 dBm. The measured value for received power from the ACS was -34.8 dBm. The 
expected (calculated) value for received power from the 12” Sphere was -47.9 dBm. The 
measured value for received power from the 12” Sphere was -46.7 dBm. All expected 
(calculated) values very nearly equaled the measured values. 

Figure 11 is from the data in the frequency domain of both the sphere and 
calibrator previously presented, now used to calculate the equivalent RCS of the 
calibrator. Measurements were made from 9-10 GHz. The VNA recorded the magnitude 
of the returned amplitude power in logarithmic form (log MAG) in units of milliwatts 
(dBm). All measurements are in the frequency domain gated for time (FG). This 
includes both the sphere (SFG) and calibrator (CFG). As can be seen, the measured 
values were very close to the expected (calculated) values. The received power for the 
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12” Aluminum Sphere was calculated to be -47.87 dBm and observed to be near -47 
dBm. The received power for the calibrator was calculated to be -36.50 dBm and 
observed to be near -35 dBm. The RCS of the sphere was previously calculated to be 
-11.37 dBsm (S-RCS). To calculate the RCS of the calibrator (C-RCS) based on the 
observed/measured data, the CFG - SFG=A (dB) was added to the sphere’s RCS. Over 
the entire frequency band the C-RCS varied from 0.6 to 3.2 dBsm. Note all these 
measurements in the chamber included the use of the RF amplifier within the calibrator 
system. When used with the weather radar, signal levels are great enough and ranges 
short enough such that an RF amplifier is not needed. Furthermore, calculations are 
made to make sure the signal return to the radar is in the proper range to prevent receiver 
saturation. 


Equivalent RCS of Calibrator 
(with in-line amplifier) 



Frequency (GHz) 


Figure 11. Frequency vs. Equivalent RCS of the ACS. (RCS_plot) 

With these measured values the equivalent RCS of the ACS was observed 
to range from 0.6 to 3.2 dBsm across the 9-10 GHz RF band. On average (mean), the 
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equivalent RCS of the ACS was 1.4934 dBsm. The median was 1.5101 dBsm. The high 
value was 3.1742 dBsm at 9 GHz and the low value was 0.5862 dBsm at 9.38 GHz. 

In summary, the basic radar measurement is echo signal power (for each 
signal sample). From the samples, average power is estimated. Reflectivity follows 
directly from substitution of average power in the weather radar equation. The ACS 
provides a method of collectively determining the imprecisely known radar system 
parameters appearing in the weather radar equation. Amplitude statistics from the 
Weather Signal Simulator, like Doppler statistics, are a planned addition once the 
programmable attenuator arrives and is installed. See Chapter IV on Future Work for 
other planned additions. 
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IV. CONCLUSION AND FUTURE WORK 

A. CONCLUSION 

An active calibration system (ACS) for weather radar was successfully designed 
and implemented. The ACS was verified to be a functional weather radar active 
calibration system under laboratory testing. It successfully provided the means (met the 
requirements) to verify all three of the signal parameters measured (estimated) by a 
weather radar. These include the zeroth Doppler moment (reflectivity), first Doppler 
moment (average radial velocity), and second Doppler moment (velocity spread). 

The cross-section (a, or equivalent RCS) of the ACS and therefore the means to 
measure and calculate the radar calibration constant (for reflectivity according to 
Equation (2.31)) was successfully determined. The underlying data to calculate the 
cross-section was measured and based on subsequent calculations determined to vary 
from 0.6 to 3.2 dBsm across the frequency band of interest. 

A constant Doppler shift (radial velocity) created by the ACS was observed and 
verified by the visual inspection of several independent human operators using a 
spectrum analyzer. The ACS successfully produced a known, operator selected radial 
velocity. While it was not possible to observe a variable Doppler shift (velocity spread) 
in exactly the same manner because of time scale limitations, the method to produce a 
variable Doppler shift was manually verified in a similar manner. The ACS produced an 
observable constant Doppler shift when provided a variable Doppler data file containing 
several entries of the same Doppler shift rather than several entries each of a different 
Doppler shift. If the control software read and executed each Doppler shift entry 
(regardless of whether they were all the same) and produced the corresponding Doppler 
shift (which could be observed), there is no reason to believe the control software would 
not read and execute each Doppler shift entry (when they are different) to potentially 
produce the corresponding varying Doppler shifts that are desired in the same manner as 
was demonstrated for the test fixture. Based on this, the ACS was deemed to have 


successfully produced a known velocity spread. 
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When tested with the MWR-05XP weather radar this should be suffieient to 
remove system bias from the refleetivity estimate and to eheek the aeeuraey of mean 
radial velocity and velocity spread estimates. 

B. FUTURE WORK 

The first item for future work is to field test the active calibrator system (ACS) 
with the MWR-05XP. It should be tested to see if the constant Doppler generated by the 
ACS is observed by the MWR-05XP. Then, a variable Doppler should be generated by 
the ACS and to see if the MWR-05XP observes a Doppler spread. It should also be 
tested to see if the equivalent RCS/Refiectivity measured is observed as expected and 
then compared to the radar calibration constants obtained from the balloon and sphere 
calibration system. 

The active calibrator system was not entirely fabricated when this thesis research 
was concluded. The programmable attenuator had not yet been received. Once received, 
the programmable attenuator will need to be installed. New RF cables will have to be 
fabricated to route the RF path sequentially through the phase shifter and programmable 
attenuator. A control/power cable will have to be fabricated to link the CRIO’s DIO 
module and the DC power supply voltages from the voltage regulator circuit with the 
programmable attenuator. This cable will be the same type (25-pin to 15-pin) as the one 
used for the same purposes by the phase shifter differing only by pin assignments. The 
parts to fabricate this cable have already been procured and received. 

The control software for the programmable attenuator has, for the most part, 
already been written. The control software to independently control the programmable 
attenuator is the same as for the phase shifter. It need only be duplicated with minor 
modifications (e.g. filename for Reflectivity values vs. filename for Doppler values, 
MATLAB script to convert Reflectivity to amplitude shifts vs. Doppler values to phase 
shifts, etc.). All of the Lab VIEW code could be duplicated with minor modifications and 
run separately to control the programmable attenuator in the same manner as the phase 
shifter (both 8-bit). In its current state, the ACS can calibrate the MWR-05XP for a 
single constant amplitude return (Reflectivity) and multiple variable Doppler values. 
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With the addition of duplicate control code and the addition of the programmable 
attenuator, the ACS can be configured to calibrate the MWR-05XP for multiple variable 
amplitude returns to simulate the changes in pulse amplitude that would be observed. 
However, controlling these values independently (one and then the other) does not allow 
a calibrator to simulate “real world” weather conditions. In order to calibrate the radar 
with “real world” weather conditions amplitude and phase values must be controlled 
simultaneously and dependently. To do this, the control code for both the programmable 
attenuator and phase shifter must be timed together. This can be achieved by controlling 
the programmable attenuator and phase shifter from within the same control code. This 
is slightly less straight forward than the minor modifications for independent control of 
the two devices. However, it is possible and desirable. The existing phase shifter control 
code would have to be modified to contain duplicate control subroutines for the 
programmable attenuator from within its own control code. The difference would 
primarily be the additional use of nested loops and sequence frames from with the same 
control code (LabVIEW VI). With implementation of this simultaneous dependent 
control “real world” weather signal statistics could be simulated by the phase and 
amplitude variation of the ACS. 

Testing of these “real world” weather patterns with the MWR-05XP is the final 
objective. However, the ACS could and should be tested incrementally with the MWR- 
05XP at various times as new features are implemented. It is suggested that the MWR- 
05XP and ACS be tested together first, “as-is” (with the phase shifter installed only), then 
after the addition of the programmable attenuator and added ability to control both the 
phase shifter and programmable attenuator simultaneously. Final testing would occur 
after all features are fully implemented and integrated into the ACS (to include output of 
amplitude and Doppler statistics of “real world” weather patterns from the Weather 
Signal Simulator). 

In summary outline. Future Work should include; 

1. A test of the ACS with the MWR-05XP in its current state (at the close of 
this thesis research). 
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a. To include observations of the first and seeond Doppler moments, 
wind speed and wind shear a.k.a. radial veloeity and veloeity 
spread a.k.a. a constant Doppler and a variable Doppler. 

b. To include an observation of the zeroth Doppler moment, a.k.a. 
reflectivity (to include use of the radar calibration constant) or a 
single, eonstant amplitude return (corresponding to the equivalent 
RCS of the ACS). 

c. To include a comparison of the observation made in (b.) to the 
observation made with a balloon and sphere calibration system. 

2. Ineorporation of the phase invariant programmable attenuator, both the 
hardware and independent software control (i.e. not affecting the digital 
phase shifter’s software control). 

3. Creation of dependent software eontrol of both the programmable 
attenuator and phase shifter (i.e. control of phase and amplitude 
simultaneously). 

4. A test of the ACS with the MWR-05XP in its advanced state (at the close 
of Future Work (2.) and (3.)). 

a. To include observations of the zeroth Doppler moment, 
reflectivities (to include use of the radar calibration constants) or 
multiple, variable amplitude returns (corresponding to several 
equivalent RCSs of the ACS). 

b. To inelude observations of all three Doppler moments 
simultaneously (simultaneous phase and amplitude). 

5. Modification of the Weather Signal Simulator to interaet directly with the 
ACS control software to output “real world” weather patterns for the 
MWR-05XP (within hardware and software limitations, both for the 
MWR-05XP and ACS - e.g. ambiguous velocities due to PRF and 
switching speed of ACS hardware). 
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6. A test of the ACS with the MWR-05XP in its final state (at the elose of 
Future Work (5.)). Are “real world” weather signal statisties observed? 

In elosing, it should also be noted that if the ACS is used at ground level (the most 
likely ease), a ground reflection will occur. The signal received by the radar will then be 
determined by the amplitudes and relative phase between the direct (line-of-sight) and 
ground reflected signal components as influenced by the ground characteristics. The 
two-way propagation gain is 12 dB if the ground reflected signal is equal in amplitude 
and in phase with the direct signal. This effect has been observed in measurements made 
with another calibration system subsequent to the completion of the research work 
described here. Compensating for this effect complicates the calibration procedure and 
will require that the ACS antenna be moved through a sequence of heights above ground 
to determine locations where the direct and ground reflected signals are in phase and 180 
degrees out of phase. Signal measurements at these two locations will permit 
determination of the propagation gain for the case where the direct and ground reflected 
signals are in phase. 
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APPENDIX A. ACS PICTURES 



Figure 12. Calibrator Front View with RAM (straight on). 



Figure 13. Calibrator Front View with RAM (looking down angle). 



Figure 14. Calibrator Left Side View (straight on). 


53 























Figure 15. Calibrator Left Side View (looking down angle). 
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Figure 16. Calibrator Left Side View (looking down angle), Zoomed. 
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Figure 17. Calibrator Left Side View (looking down angle), Zoomed - Top Cover 

OFF. 
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Figure 18. Calibrator Top View (straight on), Zoomed - Top Cover OFF. 



Figure 19. Calibrator Inside View (CRIO). 
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Figure 20. Calibrator Inside View (Power Supply Circuit). 



Figure 21. Calibrator Inside View (Phase Shifter). 
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Figure 22. Calibrator Inside View (Power Supply). 
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Figure 23. Calibrator Right Side View (looking down angle) Zoomed, Top Cover 

OFF 
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Figure 24. Calibrator Right Side View (looking from above at external in-line 

amp). 



Figure 25. Calibrator Rear Left Side View (looking down angle). 
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Figure 26. Right Horn Antenna (close-in view, polarization is horizontal). 
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Figure 27. 12” Sphere (used to calibrate the ACS). 
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APPENDIX B. DIGITAL PHASE SHIETER DATA 


Digitally Controlled Phase Shifter 
Manufaeturer; G.T. Microwave, Inc. 

Model Number; P8P-68N-5 

Electrical Specifications for Digitally Controlled 360° Phase Shifters: 
Frequency Range (GHz); 6.0-18.0 
Phase Error Vs Freq. Max: +/- 10.0° 

Amplitude Balance Max: +/- 1.50 dB 
Insertion Foss Max; 12.0 dB 
V.S.W.R. Max; 1.90:1 
Switching Speed nSec Max: 500 
RF Input Power dBm CW: +15 
RF Input Power dBm Max; +20 


Power/Fogic Connections: 


No. of 

FOGIC PIN 

+I5V 

-I5V 

GND 

BITS 

ASSIGNMENTS 

PIN 

PIN 

PIN 

8 

F.S.B. @ I to 
M.S.B. @ 8 

13 

14 

15 


Table 6. Power and Fogic Connections for Digitally Controlled Phase Shifter. 


65 





Figure 28. Specifications Sheet for Digital Phase Shifter from G.T. Microwave. 


RF Measurements on Digital Phase Shifter; See Appendix on Scalar and Vector Network 
Analyzer Measurements. 
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APPENDIX C. SCALAR AND VECTOR NETWORK ANALYZER 

MEASUREMENT DATA 


1og MflG 


MARKER 1 
3.(201 GHz 
' 35.705 dS 


,3.5015 GHz 
■ 35.003 dE 


MARKER 3 
10.002 GHz 
I 35.053 dE 


►021 

REF 35.0 dB 
2 0.5 dB.^ 

V 35.003 dB 



-i (- 


START 0.000000000 GHz 
STOP 12.400000000 GHz 


! 31 JUL 071 
L11:24:27 


Figure 29. RF Amplifier Gain Single VNA Plot. 


67 













































(11111111). 


68 





























































Figure 31. Digital Phase Shifter Insertion Loss (Left) and Phase Value (Right) Dual 
VNA Plots. Zero volts applied to MSB pin only, five volts others. 

(10000000). 
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STRRT 9.000000000 GHz 
STOP 10.000000000 GHz 


27 JLX. 07| 
17:21:46 


Figure 32. Digital Phase Shifter Insertion Loss (Left) and Phase Value (Right) Dual 
VNA Plots. Zero volts applied to second MSB pin only, five volts others. 

(01000000). 
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Figure 33. Digital Phase Shifter Insertion Loss (Left) and Phase Value (Right) Dual 
VNA Plots. Zero volts applied to third MSB pin only, five volts others. 

(00100000). 
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STOP 10.000000000 GHz 


r 27 JLL'“0- 

i 17;15: 34 


Figure 34. Digital Phase Shifter Insertion Loss (Left) and Phase Value (Right) Dual 
VNA Plots. Zero volts applied to fourth MSB pin only, five volts others. 

(00010000). 
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Sgi log MOG 

REF -10.0 dB 
0.S dB/ 


REF 10.5 « 

300.0 m °/ 


phase 


4 ^ 

c 



STPlRT B.000000000 GHz 
STOP 10.000000000 GHz 


37 JUL 0' 
17:13:06 


Figure 35. Digital Phase Shifter Insertion Loss (Left) and Phase Value (Right) Dual 
VNA Plots. Zero volts applied to fifth MSB pin only, five volts others. 

(00001000). 
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Figure 36. Digital Phase Shifter Insertion Loss (Left) and Phase Value (Right) Dual 
VNA Plots. Zero volts applied to sixth MSB pin only, five volts others. 

(00000100). 


74 

















































































Sgj_ lag MhG ►Sg^^MG phase 

REF -10.0 dB REF 2.5 ° 

0.2 dB/ 100.0 m°/ 


C 





STRRT B.000000000 GH2: 
STOP 10.000000000 GHz 


r'S? XlL 07 
[__17: 03: 5a_ 


Figure 37. Digital Phase Shifter Insertion Loss (Left) and Phase Value (Right) Dual 
VNA Plots. Zero volts applied to seventh MSB pin only, five volts others. 

(00000010). 
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0.2 dB^ 
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log MOG 


►Sgi^MS 

REF 1.5 - 

100.0 m°' 


phase 


D 



STRRT 9.000000000 GHz 
STOP 10.000000000 GHz 


! 27 JUL 07 
i 1G ;55: 37 


Figure 38. Digital Phase Shifter Insertion Loss (Left) and Phase Value (Right) Dual 
VNA Plots. Zero volts applied to LSB pin only, five volts others. 

(00000001). 
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REF -10.0 dB 

0.2 dB/ 
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►Spi/MG 

REF 0.0 “ 

100.0 m°.- 


phase 



STPRT 9.000000000 GHz 
STOP 10.000000000 GHz 


27 JUL 07 
1S:5B:49 


Figure 39. Digital Phase Shifter Insertion Loss (Left) and Phase Value (Right) Dual 
VNA Plots. Five volts applied to all pins. (00000000). 
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STRRT 9.000000000 GHz 

STOP 10.000000000 GHz 

E7 JUL 07! 
14:45:58 1 


Figure 40. Digital Phase Shifter Forward (left) and Reverse (right) Insertion Loss 
Dual VNA Plots. Demonstration of deviee’s bidireetional property. Losses 
in the forward and reverse are equal. Five volts applied to all pins. 

( 00000000 ). 
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Figure 43. Digital Phase Shifter SNA Plot. Fixed Phase State, no Doppler Shift. 
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Figure 44. Digital Phase Shifter SNA Plot. Periodic Rotation through Phase States, 

Doppler Shift. 
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Figure 46. Digital Phase Shifter Phase Value Single VNA Plot. Zero volts applied 
to second MSB pin only, five volts others. (01000000) 90° Phase Shift. 
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APPENDIX D. PROGRAMMABLE ATTENUATOR DATA 


Phase Invariant Programmable Attenuator 
Manufacturer; American Microwave Corporation 
Model Number; PI-DVAN-8012-8 

Electrical Specifications; 

Frequency; 8 GHz to 12 GHz 
Range/Resolution; 60 dB/0.25 dB 
Insertion Loss; 4.0 dB Typical 
Switching speed; 300 ns Maximum 

Maximum Power (No Damage); +27 dBm (operating +10 dBm) 
Phase Variation/Shift; 0 to 20 dB ; +/-5° 

20 to 40dB ;+/-10° 

40 to 60 dB ; +/-20° 


Power/Logic Connections; 


No. of 

LOGIC PIN 

+15V 

-15V 

GND 

BITS 

ASSIGNMENTS 

PIN 

PIN 

PIN 

8 

L. S.B. @ 5 to 

M. S.B. @ 12 

13 

14 

4 


Table 7. Power and Logic Connections for Programmable Attenuator. 


85 




ccaaows 


D ESCaiPTlOK 

*MC UOOEL ft-OV*N-8032-8 S A HIGH SPEED. lOW WSEPTION LOSS. 
phase invariant. 0 to 60 dG, S 9 iT OlOiTAL VARIAei.C. SOLID STATE. 

OCTAVr ewD ATrCNUATOR that OPERATCS over the 8 GM2 to 12 GH* RANGE. 


ORIGINAL RELEASE 


• TREOUENCV;.- .8 (Wz TO 17 CHl 

• HANGE/RESOlunON;.60<lB/'fl74<JS 

0 INSERTION LOSS;- • • .. 4.0 Cl8 TYPOI.. 3.80 <)B COAt, A.S UB MAXIMUM 

« DIGITAL CONTROL:'.8 ST TTl. 

« SWfCHlNG SPEED;.300 ns UAXHAJM 

c maximum POWER:..-i-?? 0Sm . {{DERATING 4 10 ceml 

(NO OAUACE) 

' OTOI7 ■ .36 dB'n MINMUM 

• 0IAS VOlTACSE:-.•■»12V O 200 mA MAXIMUM.. '50 mA TYPICAL 

-SV O TOO mA UAAMUU., 76 mA TYPICAL 

• PHASE VMPIATION/SHfT; qj 

O TO 20 <18.4¥ (CCMt tA’) 4 

30 TO 40 0B.ilC (CCML ±8') 

40 TO €0 40.iTO’ (GOAL tIS’l -j— 

4 Tew>ewruRe.. .q-c to tot 

• RE CONflECTORS:.SWA f£MAUE 

• BIAS VOLTAGE 

AND CONTROL. .13 PK MULTlPIH CONNECTOR 

• SI2C: ..-2.00 (L) * 2.60 (W) x O.SO (H) 2 

« WEIGHT..-< 4 OUNCES 



13 PIN MICRO D 
MULTI PIN CONtCCTOR 


40080 THRU 
4 PLACES- 



SOL® STATE 
VARIABLE ATTENUATOR 
JI unon MO. fl-PV|W-e0l2-e 



(M5 "l 

0,095 — U- 


' SNA CF> REMOVABLE 
-RF CONNECTOR 
2 PLACES 


CONFIDENTIAL AND PROPRIETARY 


ENVIRONMENTAL RATINGS : 


, ...__ <rc TO 1-70 C (OPCRATiNC) 

-65-C TO -frUSC <STORA&) 

• HCMiOnY< - . , . ... Mit-STO-20aF, MOHOO IQ3a CONO. 8 

• SHOCKS . . , . .MtL-SI0--2O2F, MCTmOO 2136 CONO. 0 

» VlGflATlCNi .... -.MIL-ST0-202E. METHOD 7040 CONO. 8 

» ALTITUDE:.MIL-ST0-202P. METHOD I05C COND. 8 

• TEMPERATUftE CYCLE' - - • U:L-STO-2<12F. HCTHOO 107D CONO- A 


ALL DIMENSIONS >IPE IN INCHES 
TOLERANCES; 

X.XX :AO.Q20 
X.XXX ±0010 


f I AMERICAN MICROWAVE CORPORATION 
LcRI PRQDCRICK, MARYLAND 


“J PRODUCT fEATuRE 

PI-0VAN-S012-8 

Jo« SOLID state: phase INVA^ANT ATTENUATOR 


Figure 47. Specifications Sheet for Programmable Attenuator from American 

Microwave. 
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APPENDIX E. ANTENNA DATA 


Antenna List: 

(2) Pasternack Sectoral Horn Antennas (used as part of ACS) 

(2) Unknown X-band Quad-Ridge Dual-Pol Horn Antennas (used as part of anechoic 
chamber ACS testing) 

(1) Narda Standard Gain Horn Antenna (used as part of anechoic chamber antenna 
testing) 
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pASTtnMACK eMTEmpnoES, me. 


PO Box 16759. Irvine. Ca 92623 


Toll Free: (866) 726-8375 
Direct: +1 (949) 261-1920 
FAX: +1 (949)261-7451 
Email: techsupport@pastemack.com 


PE9850 thru PE9864 
Standard Gain Horns 


Standard Gain Horns 

P/N 

Size 

Frequency 

(GHz) 

Gain 

(dB) 

3 dB Width 
E-Plane' H-Plane 

PE9850-10 

WR28 

26.5-40.0 

10 

54.2” 

54.4° 

PE9850-15 

WR28 

26.5-40.0 

15 

32.1° 

31.3° 

PE9850-20 

WR28 

26.5-40.0 

20 

16.7° 

18.3° 

PE9851-10 

WR34 

22.0 - 33.0 

10 

54.1° 

53,2° 

PE9851-15 

WR34 

22,0-33.0 

15 

23.1° 

40.8° 

PE9851-20 

WR34 

22.0-33.0 

20 

17.0° 

17.4° 

PE9852-10 

WR42 

18.0-26.5 

10 

58.0° 

57.0° 

PE9852-15 

WR42 

18.0-26.5 

15 

31.3° 

31.5° 

PE9852-20 

WR42 

18.0-26.5 

20 

17.5° 

17.8° 

PE9853-10 

WR51 

15.0-22.0 

10 

55.1° 

54.2° 

PE9853-15 

WR51 

15.0-22.0 

15 

32.0° 

31.8° 

PE9853-20 

WR51 

15.0-22.0 

20 

16.9° 

18.0° 

PE9854-10 

WR62 

12.4-18.0 

10 

55.3° 

50.9° 

PE9854-15 

WR62 

12.4-18.0 

15 

30.1° 

31.2° 

PE9854-20 

WR62 

12.4-18.0 

20 

18.8° 

18.9° 

PE9855-10 

WR75 

10.0-15.0 

10 

50.2° 

49.2° 

PE9855-15 

WR75 

10.0-15.0 

15 

35.4° 

28.5° 

PE9855-20 

WR75 

10.0-15.0 

20 

16.3° 

17.2° 

PE9856-10 

WR90 

8.20-12.4 

10 

48.5° 

47.4° 

PE9856-15 

WR90 

8.20-12.4 

15 

29.3° 

29.0° 

PE9856-20 

WR90 

8.20-12.4 

20 

16.1° 

16.5° 

PE9857-10 

WR102 

7.00-11.0 

10 

55.5° 

54.1° 

PE9857-15 

WR102 

7.00-11.0 

15 

29.6° 

29.3° 

PE9857-20 

WR102 

7.00-11.0 

20 

17.0° 

16.7° 

PE9858-10 

WR112 

7.05-10.0 

10 

56.8° 

55.2° 

PE9858-15 

WR112 

7.05-10.0 

15 

32.4° 

32.0° 

PE9858-20 

WR112 

7.05-10.0 

20 

19.3° 

19.3° 

PE9859-10 

WR137 

5.85-8.20 

10 

55.1° 

54.2° 

PE9859-15 

WR137 

5.85 - 8.20 

15 

33.7° 

33.2° 

PE9859-20 

WR137 

5.85 - 8.20 

20 

18.7° 

18.8° 

PE9860-10 

WR159 

4.09-7.05 

10 

59.8° 

48.3° 

PE9860-15 

WR159 

4.09-7.05 

15 

31.3° 

30.8° 

PE9860-20 

WR159 

4.09 - 7.05 

20 

14.3° 

16.9° 

PE9861-10 

WR187 

3.95-5.85 

10 

55.0° 

54.1° 

PE9861-15 

WR187 

3.95-5.85 

15 

33.8° 

33.3° 

PE9861-20 

VVR187 

3.95-5.85 

20 

18.9° 

19.2° 

PE9862-10 

WR229 

3.30 - 4.90 

10 

58.6° 

51.9° 

PE9862-15 

WR229 

3.30 - 4.90 

15 

32.7° 

32.4° 

PE9862-20 

WR229 

3.30-4.90 

20 

17.1° 

16.7° 

PE9863-10 

WR284 

2.60-3.95 

10 

50.8° 

54.1° 

PE9863-15 

WR284 

2.60 - 3.95 

15 

31.0° 

30.6° 

PE9863-20 

WR284 

2.60 - 3.95 

20 

17.2° 

16.5° 

PE9864-10 

WR430 

1.70-2.60 

10 

64.8° 

45.4° 

PE9864-15 

WR430 

1.70-Z60 

15 

33.1° 

32.0° 

PE9864-20 

WR430 

1.70-2.60 

20 

17.3° 

17.4° 



- lypicsl E-Plane Pattern 

-Typical H-Plane Pattern 

Typical 10dB Gain 



typical E-Plane Pattern 
-typical H-Plane Pattern 


Typical 15dB Gain 



-Typical E-Plane Pattern 

■ " — Typical H-Plane Pattern 


Typical 20dB Gain 


Gain and beamwidth data are typicals. 

Additional data, such as calibration data, is provided for an additional fee. 


Figure 48. Pasternack Sectoral Horn Antenna Data Sheet. 
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Figure 49. Pasternack Sectoral Horn Antenna Outline Drawing. 
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NOTES: 

1. UNLESS OTHERWISE SPECIFIED ALL DIMENSIONS ARE NOMINAL. 

2. ALL SPECIFICATIONS ARE SUBJECT TO CHANGE WITHOUT NOTICE AT ANY TIME. 





2.60-40 GHz 

Standard Gain Horns 

• Primary Standard of Antenna Gain 
^ j lyiQjjgIg 2.60 GHz to 40 GHz 



FREQUENCY 

RANGE 

(GHz) 

MODEL 

SAND 

WAVEGUIDE 

SIZE 

INPUT 

COVER FLANGE 
EQUIVALENT 

VSWR 

(Max) 

WEIGHT 

Iba, Kg 

2.60-3.95 

644 

S 

WR-284 

UG-584/U 

1.15 

6 

2,8 

: .-h. .. . 

643 

c 

WR-187 

UG-407/U 

1.15 

2.3 

1.1 

5.4-8.20 

642 

XN 

WR-137 

UG-441/U 

1.15 

1.0 

0.5 

8.20-12.4 

640 

X 

WR-90 

UG-135/U 

1.15 

0.5 

0,23 

12.4-18 

639 

KU 

WR-62 

UG-419/U 

1.15 

0.2 

0.1 

18-26.S 

638 


WR-42 

UG-595/U 

1.15 

0.2 

0,1 : 

26.5-40 

V637 

V 

WR-28 

UG-599/U 

1.15 

0.1 

0.05 


Patterns for al models in this series conform to the following 
description; Beam width in E and H plane varies from 23'' at the 
highest frequency to SA-'-' at the lowest frequency. Side lobes in 
the H plane are all more than 20 dB down. First side oPes in the 
E p ane are 13 dB down, second side lobes are 18 dB down and 
all other E plane iooes are more than 20 dB down. 

Gain at Mid Frequency; 16.5 dB (with reference to isotropic radiation) 
variation is 1.5 dB over tota' band about the mid band va ue. 


narda 


435 MORELAND ROAD - HAUPPAUQE. NY 11788 
USA "EL: (1) 631 231-1700 • INT'LTEL: (1) 63 1 231-1390 • FAXri) 631 231-1711 
E-MAIL: naraaeasT©L-3COM,Gom * wvvw.narGamicrowave.com 


Figure 51. Narda Standard Gain Horn Antenna Datasheet (page 1). 
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Waveguide Horns 


Standard Gain Horn 






.i-A.i-i. 

. [-■r-B-- . 

': 'MAX:'' 7 '':' 

.. 

1 'MAX :. 

V637 

1.76 

1.06 

.82 

(44,7) 

(26.9) 

(20.8) 

636 

2.57 

1.51 

1.16 


(36.3) 

(29.4) 

639 

3.48 

2.20 

1.73 

/3) 

(55.8) 

(43.9) 

640 

5.06 

3.09 

3.34 

(128.5) 

(78.4) 

(53,4) 

642 

7,76 

4,67 

3.53 

(197.1) 

(118.6) 

(83,6) 

B43 

10.47 

6.34 

4.80 

(265.9) 

(161.0) 

(121.9) 

644 

.62 

9.52 

7.16 

(401.8) 

(241.8) 

(181,8) 



Dimensions in oarentheses are in miilimeters and tor reference ony. 
■ ■: ■ A'l dimensions max, un;ess othem/ise specified/ ■: 


435 MORELAND ROAD • HAUPPAUQE. NY 1178S 

USA TEL: (I) 631 231-1700 ® !N~LTEL: (1) 631 231-1330 • FAX(I) 631 231-1711 
E-MAIL: nRrdReafit©L-3COM,Qorr * wvvv.'.narnaiTiicroy.'ave.com 


narda 



Figure 52. Narda Standard Gain Horn Antenna Datasheet (page 2) 














































WAVEGUIDE TO COAXIAL ADAPTERS FOR STANDARD GAIN HORN 


STANDARD GAIN 

HORN MODEL 

644 

; 643 : 

WAVEGUIDE TO COAXIAL TYPE 

N - FEMALE MODEL 

614A 

613A ; : 

ADAPTER SIWA/3.5/2.9 
FEMALE MODEL 

642 

612A 

4602 




.. 

■ : ■ : 60TA: ■ 

: : : 4601] : 1 : 1 : 

639 

609 

4609 







. . .. . . . 


. . .. . 

V637 


V4607 



narda 


435 MORELAND ROAD - HAUPPAUQE. NY 11788 
USA "EL: (1) 631 231-1700 • INT'LTEL: (1) 63 1 231-1390 • FAXri) 631 231-1711 
E-MAIL: naraaeasT©L-3COM,Gom * www.naraamicrowave.coifi 


Figure 53. Narda Standard Gain Horn Antenna Datasheet (page 3). 
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ABSOLUTE OAiN OALieBATIQN 

H0OiL S4S STANDARD 6AIN H0RN 



t.f i.p i».e 11.9 11.9 i:.e is.? 

fHEGUEHCT fOnZI 


AESOLUTS OAEN OALiSRAtlON 
NARDA MODEL €43 STANDARD QAtN HORN 



i^nEsusHCT laos 


ABSOLUTS SA3M CAUBBATE9R 
NAN DA MODEL @4S STANDARD OAIN HORN 



pRsanEfseY faHsi 


narda 


435 MORELAND ROAD - HAUPPAUQE. NY 11788 
USA "EL: (1) 631 231-1700 • INT'LTEL: (1) 63 1 231-1390 • FAXri) 631 231-1711 
E-MAIL: naraaeasT@L-3COM.Gom * wv'w'w.naraamicrowave.com 


Figure 54. Narda Standard Gain Horn Antenna Datasheet (page 5). 
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APPENDIX F. ANTENNA PATTERN DATA 


A total of five antennas were used in testing. One Narda standard gain horn 
antenna was used during antenna testing (patterns) only. Two X-band horn antennas 
were used during antenna testing (patterns), then used during ACS testing (as the transmit 
and receive antennas of the radarA^NA). Two Pasternack horn antennas were used 
during antenna testing (patterns), then used during ACS testing (as the transmit and 
receive antennas of the ACS). Measurements were made at three frequencies: 9370, 
9680, and 9990 MHz (low, mid, and high frequency points on the MWR-05XP weather 
radar’s RF band). E-plane and H-plane patterns were measured. To measure both 
planes, the antenna under test was first rotated 360° in azimuth with the radiator oriented 
for horizontal polarization (E-plane), then rolled 90° and the process repeated with the 
radiator now oriented for vertical polarization (H-plane). In total, two plane patterns 
were measured for five antennas at three frequencies to create thirty data sets. 

A. VNA ANTENNA PATTERN DATA 

In all plots, solid line is for X-band Horn Antenna #1 (“xbh”) and dotted line is 
for X-band Horn Antenna #2 (“ybh”). The last three plots are solid line for the Standard 
Gain Horn (“sgh”) which are included to provide a gain reference. 
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Figure 55. X-band Horn Antenna #1 Pattern at 9370MHz (E plane in red and H 

plane in blue). 



Figure 56. X-band Horn Antenna #2 Pattern at 9370MHz (E plane in red).i2 


12 Note because the measured values slightly exceeded -15 dB the gain values have been normalized to 
0 rather than -15. 
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Figure 57. X-band Horn Antenna #2 Pattern at 9370MHz (H plane in blue). 



Figure 58. X-band Horn Antenna #1 Pattern at 9680MHz (E plane in red and H 

plane in blue). 
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Figure 59. X-band Horn Antenna #2 Pattern at 9370MHz (E plane in red and H 

plane in blue). 



Figure 60. X-band Horn Antenna #1 Pattern at 9990MHz (E plane in red and H 

plane in blue). 
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Figure 61. X-band Horn Antenna #2 Pattern at 9990MHz (E plane in red and H 

plane in blue). 



Figure 62. Standard Gain Horn Antenna Pattern at 9370MHz (E plane in red and H 

plane in blue). 
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Figure 63. Standard Gain Horn Antenna Pattern at 9680MHz (E plane in red and H 

plane in blue). 



Figure 64. Standard Gain Horn Antenna Pattern at 9370MHz (E plane in red and H 

plane in blue). 
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B. ACS ANTENNA PATTERN DATA 


In all plots, solid line is for Pasternack Horn Antenna #1 (a.k.a. “xpe”) and dotted 
line is for Pasternack Horn Antenna #2 (a.k.a. “ype”). The last three plots are solid line 
for the Standard Gain Horn (a.k.a. “sgh”). The pattern plots for the Standard Gain Horn 
have been included to provide a gain reference and in all three plots the gain has been 
adjusted to account for the use of a 6.5 dB attenuator. Also, note all data have been 
shifted -90° to boresight the patterns at 0°. This was done to make these patterns look 
like the other patterns previously presented. Additionally, in all plots the E-plane pattern 
is in red and H-plane pattern is in blue. 



Figure 65. Pasternack Horn Antenna #1 Pattern at 9370MHz (E plane in red and H 

plane in blue). 
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Figure 66. Pasternack Horn Antenna #2 Pattern at 9370MHz (E plane in red and H 

plane in blue). 



Figure 67. Pasternack Horn Antenna #1 Pattern at 9680MHz (E plane in red and H 

plane in blue). 
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Figure 68. Pasternack Horn Antenna #2 Pattern at 9680MHz (E plane in red and H 

plane in blue). 



Figure 69. Pasternack Horn Antenna #1 Pattern at 9990MHz (E plane in red and H 

plane in blue). 
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Figure 70. Pasternack Horn Antenna #2 Pattern at 9990MHz (E plane in red and H 

plane in blue). 



Figure 71. Standard Gain Horn Antenna Pattern at 9370MHz (E plane in red and H 

plane in blue). 
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Figure 72. Standard Gain Florn Antenna Pattern at 9680MHz (E plane in red and H 

plane in blue). 



Figure 73. Standard Gain Horn Antenna Pattern at 9990MHz (E plane in red and H 

plane in blue). 
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C. MATLAB PLOTTING FILES 

plotpatterns.m _ 


polarDB(xbh9370H(;, 1 ),xbh9370H(: ,2),’r') 
polarDB(xbh9370V(;, 1 ),xbh9370V(: ,2),’b’) 

figure 

%% normalized ro to 0 not -15 

polarDB(ybh9370H(;,l),ybh9370H(:,2),’r;') 

figure 

polarDB(ybh9370V(;,l),ybh9370V(;,2),'b;') 

figure 

polarDB(xbh9680H(;,l),xbh9680H(;,2),V) 

polarDB(xbh9680V(;,l),xbh9680V(;,2),’b') 

figure 

%% both normalized to 0 

polarDB(ybh9680H(:,l),ybh9680H(:,2),'r;') 

polarDB(ybh9680V(;,l),ybh9680V(:,2),’b;’) 

figure 

polarDB(xbh9990H(;, 1 ),xbh9990H(; ,2),’r') 
polarDB(xbh9990V(:, 1 ),xbh9990V(; ,2),’b') 

figure 

polarDB(ybh9990H(;,l),ybh9990H(:,2),'r:') 
polarDB(ybh9990V(;, 1 ),ybh9990V(; ,2),'b;') 

figure 

polarDB(sgh9370H(;,l),sgh9370H(;,2),'r-') 
polarDB(sgh9370V(;, 1 Xsgh9370V(; ,2),’b-') 

figure 

polarDB(sgh9680H(;,l),sgh9680H(;,2),'r-') 

polarDB(sgh9680V(;,l),sgh9680V(;,2),'b-') 

figure 

polarDB(sgh9990H(;, 1 ),sgh9990H(; ,2),'r-') 
polarDB(sgh9990V(:, 1 ),sgh9990V(: ,2),’b-') 
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plot pattems2.m 


%% -90 correction for 90 degree offset when patterns taken (-90 moves 
%% boresight to 0 degrees) 
polarDB(xpe9370H(;,l)-90,xpe9370H(:,2),'r') 
polarDB(xpe9370V(;,l)-90,xpe9370V(;,2),'b') 

figure 

polarDB(ype9370H(:,l)-90,ype9370H(;,2),'r;') 

polarDB(ype9370V(;,l)-90,ype9370V(;,2),'b;') 

figure 

polarDB(xpe9680H(;,l)-90,xpe9680H(;,2),'r') 

polarDB(xpe9680V(;,l)-90,xpe9680V(;,2),'b') 

figure 

polarDB(ype9680H(;,l)-90,ype9680H(;,2),'r;') 

polarDB(ype9680V(:,l)-90,ype9680V(:,2),'b;') 

figure 

polarDB(xpe9990H(:,l)-90,xpe9990H(;,2),'r') 

polarDB(xpe9990V(;,l)-90,xpe9990V(;,2),'b') 

figure 

polarDB(ype9990H(;,l)-90,ype9990H(;,2),'r;') 

polarDB(ype9990V(:,l)-90,ype9990V(:,2),'b;') 

%%+6.5 correction for 6.5dB attenuator used in sgh measurement, not used 

%%in pe measurements 

figure 

polarDB(sgh9370H(:,l),sgh9370H(:,2)-l-6.5,'r-') 

polarDB(sgh9370V(:,lXsgh9370V(;,2)+6.5,’b-') 

figure 

polarDB(sgh9680H(:,l),sgh9680H(;,2)-l-6.5,'r-') 

polarDB(sgh9680V(;,l),sgh9680V(;,2)+6.5,'b-') 

figure 

polarDB(sgh9990H(;,l),sgh9990H(;,2)-l-6.5,'r-') 

polarDB(sgh9990V(:,l),sgh9990V(:,2)+6.5,’b-') _ 
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polarDB.m 


function polardb(ang,pat,line_style) 

% NON-NORMALIZING VERSION THAT DOES NOT REQUIRE POEAR2 
% take angle and dB data and generate polar plot 
% angle data in ang(;) is degrees; dB data in pat(;) 
if any(size(ang) ~= size(pat)) 

error('THETA and RHO must be the same size.'); 
end 

pmax=max(pat); 

dynr=60; % dynamic range of dB scale 
rticks=4; % number of rings in polar plot 

dbstep=dynr/rticks; % dB spacing between rings — must be integer 
top=floor(pmax/dbstep)*dbstep+dbstep; % convenient outer ring dB value 
theta=ang*pi/180; 

% set a floor 
for i=l ;length(pat) 

if pat(i) < top-dynr, pat(i)=top-dynr; end 
end 

rho=p at+dynr-top; 
hold on 

% define a circle 

th = 0;pi/50;2*pi; 
xunit = cos(th); 
yunit = sin(th); 

% now really force points on x/y axes to lie on them exactly 
inds = [l;(length(th)-l)/4;length(th)]; 
xunits(inds(2;2;4)) = zeros(2,l); 
yunits(inds(l;2;5)) = zeros(3,l); 
rmax=dynr; rmin=0; 
rinc = (rmax-rmin)/rticks; 
for i=(rmin+rinc);rinc;rmax 

plot(xunit*i,yunit*i,'-','color','k','linewidth',l); 

text(0,i+rinc/10,[' ' num2str(i-dynr+top)],'verticalalignment','bottom' ); 
end 

% plot spokes 

th = (l;6)*2*pi/12; 

cst = cos(th); snt = sin(th); 

cs = [-cst; cst]; 

sn = [-snt; snt]; 

plot(rmax*cs,rmax*sn,'-','color','k','linewidth',l); 

% annotate spokes in degrees 
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rt = 1.1 *rmax; 

for i = 1 ;max(size(th)) 

text(rt*cst(i),rt*snt(i),mt2str(i*30),'horizontalalignment','center'); 
if i == max(size(th)) 

loc = int2str(0); 

else 

loc = int2str(l 80+1*30); 
end 

text(-rt*cst(i),-rt*snt(i),loc,'horizontalalignmenf,'center' ); 
end 

% set viewto 2-D 
view(0,90); 

% set axis limits 

axis(rmax*[-l 1-1.1 1.1]); 

% transform data to Cartesian coordinates 
XX = rho.*cos(theta); 
yy = rho.*sin(theta); 

% plot data on top of grid 

q = plot(xx,yy,line_style); 
hpol = q; 

axis('equar);axis('off); 
hold off 
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D. ABSOLUTE GAIN AND 3 DB BEAMWIDTH FOR X-BAND AND 
PASTERNACK HORN ANTENNAS (CALCULATIONS AND 
MEASUREMENTS) 


Standard Gain Horn Data from Manufacturer's Chart 

Frequency (GHz) Gain (dB) 

9.37 15.66 

9.68 15.90 

9.99 16.15 

Standard Gain Horn (sgh) Data Measured in Chamber 


Frequency (GHz) 

Plane 

Location (Degrees 

Gain (dB) 









(E or H) 

from Boresight) 









9.37 

E 

-1 

-12.43 








9.37 

H 

2 

-12.61 








9.68 

E 

0 

-12.07 








9.68 

H 

2 

-12.65 








9.99 

E 

0 

-12.05 








9.99 

‘slight variation in setup 

H 

0* 

-11.75 








X-band Horn 1 (xbh) 

M=Measured, C=Calculated 

M 

C 

C 

M 

M 

M 

M 

C 

Frequency (GHz) 

Plane 

Location (Degrees 

Relative Gain 

A sgh (dB) 

Absolute 

First 

-3dB Point 

Second 

-3dB Point 

3dB Beamwidth 

(Eor H) 

from Boresight) 

(MAX Value, dB) 

Gain (dB) (Loc, Deg) 

(Rel Gain, dB) 

(Loc, Deq) 

(Rel Gain, dB) 

(Degrees) 

9.37 

E 

2 

-15.18 

2.75 

12.91 

-15 

-18.04 

17 

-18.08 

32 

9.37 

H 

3 

-15.66 

3.04 

12.61 

-21 

-18.64 

24 

-18.62 

45 

9.68 

E 

3 

-15.10 

3.02 

12.88 

-14 

-18.09 

16 

-18.03 

30 

9.68 

H 

2 

-15.11 

2.46 

13.44 

-19 

-18.14 

22 

-18.10 

41 

9.99 

E 

0&1 

-15.60 

3.55 

12.60 

-14 

-18.42 

16 

-18.56 

30 

9.99 

H 

1 

-15.72 

3.97 

12.18 

-18 

-18.62 

22 

-18.73 

40 

X-Band Horn 2 (ybh) 

9.37 

E 

2 

-14.71 

2.28 

13.38 

-14 

-17.43 

17 

-17.58 

31 

9.37 

H 

3 

-15.02 

2.41 

13.25 

-21 

-18.02 

24 

-18.06 

45 

9.68 

E 

3 

-14.89 

2.82 

13.08 

-13 

-17.81 

17 

-18.04 

30 

9.68 

H 

2 

-14.80 

2.15 

13.75 

-19 

-17.93 

22 

-17.74 

41 

9.99 

E 

1 

-15.12 

3.06 

13.09 

-14 

-18.09 

16 

-18.11 

30 

9.99 

H 

2 

-15.34 

3.59 

12.56 

-18 

-18.38 

22 

-18.29 

40 


Table 8. Datafile for X-band Horn Antennas. 
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Figure 74. Plot of Absolute Gain for X-band Flom Antennas. 
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Figure 75. Plot of 3dB Beamwidth for X-band Florn Antennas. 


Standard Gain Horn 

Data from Manufacturer’s Chart 








Frequency (GHz) 

Gain (dB) 










9.37 

15.66 










9.68 

15.90 










9.99 

16.15 










Standard Gain Horn (sgh) 

Data Measured in Chamber [minus -6.5dB attenuator in line] 






Frequency (GHz) 

Plane 

Location (Degrees 

Gain (dB) 









(Eor H) 

from Boresight) 









9.37 

E 

-1 

-5.93 

-12.43 







9.37 

H 

2 

-6.11 

-12.61 







9.68 

E 

0 

-5.57 

-12.07 







9.68 

H 

2 

-6.15 

-12.65 







9.99 

E 

0 

-5.55 

-12.05 







9.99 

‘slight variation in setup 

H 

0* 

-5.25 

-11.75 







Pasternack Horn 1 (1PE) 

M=Measured, C=Calculated 

M 

C 

C 

M 

M 

M 

M 

C 

Frequency (GHz) 

Plane 

Location (Degrees 

Relative Gain 

A sah fdBl 

Absolute 

First 

3dB Point 

Seconc 

-3dB Point 

3dB Beamwidth 


(Eor H) 

from Boresight) 

(MAX Value, dB) 


Gain (dB) (Loc, Deg) 

(Rel Gain, dB) 

(Loc, Deq) 

(Rel Gain, dB) 

(Degrees) 

9.37 

E 

92 

-7.24 

1.31 

14.35 

76 

-10.30 

109 

-10.32 

33 

9.37 

H 

93 

-8.85 

2.74 

12.91 

77 

-12.04 

108 

-12.10 

31 

9.68 

E 

92 

-6.87 

1.30 

14.60 

77 

-9.93 

108 

-9.91 

31 

9.68 

H 

92 

-8.16 

2.01 

13.89 

79 

-10.96 

106 

-10.96 

27 

9.99 

E 

91 

-7.02 

1.47 

14.68 

76 

-10.02 

105 

-10.19 

29 

9.99 

H 

93 

-8.43 

3.18 

12.97 

78 

-11.61 

107 

-11.57 

29 

Pasternack Horn 2 (2PE) 

9.37 

E 

93 

-7.38 

1.45 

14.20 

77 

-10.21 

109 

-10.30 

32 

9.37 

H 

93 

-8.86 

2.74 

12.91 

78 

-11.67 

107 

-11.64 

29 

9.68 

E 

93 

-6.71 

1.14 

14.76 

77 

-9.87 

109 

-9.92 

32 

9.68 

H 

92&93 

-8.07 

1.92 

13.98 

79 

-10.91 

107 

-11.24 

28 

9.99 

E 

93 

-6.85 

1.29 

14.86 

78 

-9.98 

108 

-10.19 

30 

9.99 

H 

93 

-8.17 

2.92 

13.23 

79 

-10.96 

107 

-11.27 

28 


Table 9. Datafile for Pasternack Flom Antennas. 
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Figure 76. Plot of Absolute Gain for Pasternaek Florn Antennas. 
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3dB Beamwidth 



—♦—Pasternack Horn 1 (1PE)- H-Pol/E-Field 
-■-Pasternack Horn 1 {1PE)- V-Pol/H-Field 
-♦- Pasternack Horn 2 (2PE) - H-Pol/E-Field 
Pasternack Horn 2 (2PE) - V-Pol/H-Field 


Frequency (GHz) 


Figure 77. Plot of 3dB Beamwidth for Pasternaek Florn Antennas. 


E. EXPLANATION OF DATAFILE MEASUREMENTS AND 
CALCULATIONS 

The first entry in both datafiles is the absolute gain of the referenee horn antenna, 
the Narda standard gain horn antenna. These values were obtained from the Narda 
datasheet, page 5. See below for a repliea of this graph with the look-up proeess 
annotated. 
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Figure 78. Narda Standard Gain Horn Absolute Gain Plot (zoomed to area of 
interest and annotated with values of interest). 


The second entry in both datafiles is the measured power received when using the 
standard gain horn antenna during antenna testing as the antenna under test. In the 
Pasternack datafile these values have been adjusted by 6.5 dB because the measured 
values included an additional 6.5 dB of loss from an in-line attenuator that was required 
to prevent the power level from saturating the measurement device (Vector Network 
Analyzer or VNA). 

The third entry in both datafiles is the measured power received when using either 
the X-band or Pasternack horn antenna during antenna testing as the antenna under test. 
The maximum relative amplitude as measured is listed. Then this value is subtracted 
from the corresponding standard gain horn measurement value to get a delta. Note the 
value is subtracted because the power values are in negative dBs and the X-band and 
Pasternack horn values were less than the standard gain horn values. This delta is then 
subtracted from the value of the absolute gain of the reference horn antenna obtained 
from the graph. Note the values are subtracted because the delta though positive is the 
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amount below in power. The result is the absolute gain for that antenna at the given 
frequency. 

Then the value for the first measurement to exceed 3 dB below the maximum 
relative amplitude is listed along with its angle in degrees. This was done for both the 
first value above and below in angle of the maximum relative amplitude angle. These 
angle values were subtracted from one another to find the 3dB Beamwidth. 

Example calculation from X-band Horn #1 at 9.37 GHz in the E plane (as 
described above); 

Absolute Gain; 

(-12.43115)-(-15.18115) = 2.75 dB 
15.655-2.75 = 12.905 dB 

3dB Beamwidth; 

17-(-15) = 32 dB 

F. ANTENNA ISOLATION TEST DATA 

At 80 cm separation from each other oriented main beam parallel to each other 
both at horizontal polarization the isolation between the antennas was found to be -42 dB. 
At cross polarization (one horizontal and one vertical) the isolation was found to be -57 
dB. At slant 45° degree polarization (parallel, then perpendicular) the isolation was 
found to be -47 dB. At twice the separation distance (160 cm), the isolation increased by 
6 dB (to -48 dB from -42 dB for co-polarization horizontal). This demonstrated the drop¬ 
off from spreading loss alone (where twice the distance reduces power by a quarter). A 
distance of 6 ft (same as lumber length) was chosen for the antenna separation and 
mounting on the ACS. This was sufficient isolation to prevent regenerative feedback. In 
addition to this minimum antenna separation distance an RE isolator was used to make 
the ACS transmission path regenerative. There is a need for an RE isolator because all of 
the RE components are bidirectional. When the RE amplifier was used (as in the 
anechoic chamber testing) the isolator was not required, as an amplifier is a uni¬ 
directional device. If the signal path through the ACS is allowed to be bi-directional, 
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then the horns beeome radiators in a two element array with large eleetrieal spaeing. This 
results in a radiation pattern with numerous narrow lobes. This eondition must be 
eliminated if the ACS radiation pattern is to be that of a single horn. 
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APPENDIX G. ANECHOIC CHAMBER 


This appendix provides information on the setup of the Naval Postgraduate 
SehooTs aneehoie ehamber on the roof of Spanagel Hall. This is the aneehoie ehamber 
owned by the Eleetrieal and Computer Engineering Department. This was the only 
aneehoie ehamber used in this thesis (to inelude Antenna Testing and 
Amplitude/Refleetivity Measurements). 
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Anechoic Chamber Calculations 


Distance between homs - 19 ft. about 18.5 ft for the big horn 


Path loss @ 5.8GHz = 10log 


10 


4^' 

^ c J 




General Friis Equation 

= P^-\-G,+G^- Lp-\- gains - losses 
In the chamber 

“ ^palh ~ Kariabteau ~ pad 

The chamber measures P^ - P^^ or the system loss. 
Solving for P^ and substituting into the equation for P^, 
after a little algebra, we find that; 



[p -P . 
y r ref 

measured value 


-G,+L^+L 


variable an 


-30dB 


MP«7300C 



Planned configuration after repairs 
March 2006 


Figure 79. Anechoic Chamber Diagram. 
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APPENDIX H. COMPACTRIO 



Figure 80. Figure of CompactRIO Software Architecture from ni.com. 


A. COMPACTRIO REAL-TIME EMBEDDED CONTROLLER (NI CRIO- 

9004, 64/512 MEMORY) 

This model embedded controller runs Lab VIEW Real-Time with 64 MB of 
nonvolatile CompactFlash storage and 512 MB of DRAM memory, has a 195 MFlz 
processor, and is equipped with a lO/lOOBaseT Ethernet port. 

Erom the CompactRIO ADVISOR on the webpage for this device, 
http://sine.ni.eom/nips/cds/view/p/lang/en/nid/I4I55 “National Instruments CompactRIO 
real-time embedded controllers offer powerful stand-alone embedded execution for 
deterministic EabVIEW Real-Time applications. The NI CompactRIO controllers are 
designed for extreme ruggedness, reliability, and low power consumption, with dual 9 to 
35 VDC supply inputs that deliver isolated power to the CompactRIO chassis/modules 
and a -40 to 70 °C operating temperature range. The NI cRIO-90Ix series controllers 
feature a 400 MHz industrial processor balancing low power consumption and powerful 
real-time floating-point signal analysis, control and logging. The cRIO-900x controllers 
feature a 200 MHz industrial processor for real-time analysis, control and logging.” 
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CompactRIO Real-Time Embedded Controllers 
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Owerwiew and Applications 

National Instruments cRIO-900x real-time embedded controllers offer 
powerful stand-alone embedded execution for deterministic LabVIEVV 
Real-Time applications. The Nl cRiO-9002 includes 32 MB of DRAM 
memory and C4 MB of nonvolatile flash storage for file storage, [he 
cRIO-9004 it^cluiles 64 MB of DRAM memory arid 51 ? MB of ncnvulalile 
Hash storage lor data-logging applications. Rolh contreilers arc designed 
for extreme ruggedness, reliability, and low power consumption with 
dual Oto dBVDC supply inputs that deliver isolated power to the 
CompactRIO chassis/mcdules and a -dOtc 70 '"C temperature range, 

A195 MHz industrial processor balances low pcwiTr consumption with 
powertul reaMimo lloaling-point signal processing and analysis 
capabilihes for deterministic cunlroi loops exceeding 1 kHz, 



figure 1. CGnitj<ict.RiO Harriwsre Architecture 


System Configuration 

Iho CompaclRIO real time controller conncclsto any 4 or 8 slot 
CornpaclRiO reconfiyurabie chassis. Tlie user-defined FPGA circuitry in 
the chassis ccnlrols each I/O module and passes dala lo the ccntrollor 
through a local PCI bus, using built-in communication functions. 


Embedded Software 

You can synchronize embedded code execution to an FPGA-generated 
interrupt request (IRQ) or an internal millisecond real-time clock source. 
The LabVIFW Real-Time F.TS OS provides reliatiilily and simplifies tlie 
development ol complete embedded applical.Icns [hat include l.ime-icritical 


ccnlrul and acquisition loofis in additiun to lower-priority loops for 
postprocessing, data logging, and Ethernet/serial ccmmunication. 

Built-in elemental I/O functions such as the FPGy\ Read,/Write function 
provide a ccmmunication interface to the highly optimized reconfigurabSe 
PHOA circuitry. Data values are read trem the PPb'A in integer icrmal, and 
llinii corsverled lo scaled Kigineeriny units in the ccntruiier. 
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Figure 82. Datasheet for CompaetRIO Real-Time Embedded Controllers from 

ni.com. 
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B. COMPACTRIO RECONFIGURABLE EMBEDDED CHASSIS (NI CRIO- 

9101, IM GATE 4-SLOT CHASSIS) 

From the CompactRIO ADVISOR on the webpage for this device, 
http://sine.ni.eom/nips/cds/view/p/lang/en/nid/14154 “The National Instruments 
CompactRIO embedded chassis are the heart of the CompactRIO system because they 
contain the reconfigurable I/O (RIO) FPGA core. The RIO core has an individual 
connection to each I/O module and is programmed using easy to use elemental I/O 
functions to read or write signal information from each module. Because there is no 
shared communication bus between the RIO FPGA core and the FO modules, 
input/output operations on each module can be precisely synchronized with 25 ns 
resolution. The RIO core can perform local integer-based signal processing and decision 
making and directly pass signals from one module to another.” 



Figure 83. Picture of cRIO-9103 from ni.com. 
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CompactRIO 
Reconfigurable Chassis 



• Design hardware using Lab VIEW 
programming skills 

• 4 or 8-slot chassis for any 
CompactRIO I/O modules 

• 1 M or 3 M gate RIO FPGA core 
for normal or extended RIO 
processing power 

• DIN-raiL mounting, 19 in. rack 
mount, and panel mounting options 

• NI CompactRIO Extreme Industrial 
Certifications and Ratings'- 


• Program in easy-to-use 
LabVIEW FPGA graphical 
development environment 
to automatically synthesize 
an optimized high- 
performance electrical 
circuit implementation 

of your application 

• RIO FPGA core executes 
LabVIEW control logic at 
rates up to 40 MS/s using 
single-cycle timed loops 




Product 

■ 'cRiOiaioi' 
cRIO-9102 


' cRIO/9103' 
cRIO-9104 


FPfiA System Gates 

'liM 


RAM (KB) 

; ■ 82' ■ ■ 


DefaultTimebase(MHz} 


Maximum Power Consumption (W) 


Built-In 

Panei Mounting Hoies 


Overview and Applications 

The National Instruments CompactRIO reconfigurable chassis are 
the heart of the CompactRIO system because they contain the 
reconfigurable I/O (RIO) core. The RIO FPGA core, which has an 
indi-vidual connection to each I/O module, is programmed with easy- 
to-use elemental I/O functions to read or write signal information 
from each module. Because there is no shared communication bus 
between the RIO FPGA core and the I/O modules, I/O operations on 
each module can be precisely synchronized with 25 ns resolution. The 
RIO core can perform local integer-based signal processing and 
decision-making and directly pass signals from one module to another. 

The RIO core is also connected to the CompactRIO real-time 
controller through a local PCI bus interface. The real-time controller 
can retrieve data from any control or indicator on the front-panel of 
the RIO FPGA application throt^ an easy-to-nse FPGA Read/Write 
function. The RIO FPGA can also generate interrupt requests (IRQs) 
to s-yndironize the real-time software execution -with the RIO FPGA. 
Typically, the real-time controller is used to convert the integer based 
I/O data to scaled floating-point numbers. In addition, the real-time 
controller typically performs single-point control, waveform analysis, 
data losing, and Ethemet/serial communication. 

The reconfigurable chassis, real-time controller, and I/O modules 
combine to create a complete stand-alone embedded system. 
Application development consists of three steps: 

1. Target the reconfigurable chassis to automatically detect the 
I/O modules and develop the RIO FPGA application, 

2. Compile the RIO application to automatically synthesize 

an optimized high-performance electrical circuit implementation 
of your application, 

*See CompactRIO Overview on page 366 for details. 


3. Develop the real-time controller application to add 
floating-point control, signal processing, data logging, 
and communication 

Key Features 

• Create any local or multichassis timing, tri^ering, and 
synchronization scheme with 25 ns resolution 

• Use multiple while loops to create a parallel processing application 
for high-performance signal processing or multirate control systems 

• Bufit-in PID control fimctions for control system loop rates 
greater than 100 kHz 

• Generate waveforms or implement nonlinear look-up tables 
(LUTs) using LabVIEW FPGA express Vis 

»Integrate -widely available third-party HDL cores using the 
LabVIEW FPGA Module HDL Node 

»Enforce critical logic and interlocks in siiicon hardware drcuitry, or 
use the parallel RIO architecture to create duai, triple, or quadruple 
redundant systems 

Visit m.com/cGmp&cJjiG for example programs, application notes, 

and other developer tools. 

For ordering information, see page 381. 


.®« ® d©- <g> g C€ © 

h=?sRi;ii( I xsKjr? (pending) 

For more information, see page 379. 

BUY ONLINE ot CALL '866.' 266-9891! f: ni.coni 
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Figure 84. Datasheet for CompaetRIO Reeonfigurable Embedded Chassis from 

ni.eom. 
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C. COMPACTRIO I/O MODULES FOR DIGITAL INPUT AND 
COUNTER/TIMER (NI 9401, 8 CH 5 V/TTL HIGH SPEED 

BIDIRECTIONAL DIGITAL I/O MODULE) 


From the CompactRIO ADVISOR on the webpage for this device, 
http://sine.ni.eom/nips/cds/view/p/lang/en/nid/14168 “National Instruments CompactRIO 
digital input modules are compatible with any CompactRIO chassis and work with 5 
V/TTL, industrial and universal logic levels and signals for direct connection to a wide 
array of devices such as industrial switches, encoders and transducers. Each module 
features an integrated connector junction box with a screw terminal or industrial D-Sub 
connector for flexible and low cost signal wiring. All modules feature the NI 
CompactRIO Extreme Industrial Certifications and Ratings. Each CompactRIO EO 
module has a direct connection to the reconfigurable EPGA in CompactRIO systems. 
Using the EabVIEW EPGA module, users can configure and reconfigure all CompactRIO 
digital input modules to implement counter/timer, encoder, digital protocol, and/or static 
digital inputs.” 



Eigure 85. Picture of cRIO-9401 from ni.com. 
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Owerwiew 

High-pertorniancu digital output and switching modules tor National 
Instruments CcmpactRIQ embedded systems, R Series expansion 
chassis, and Nl CompactDAQ systems provide e.xteiided voltage ranges 
and high current-switching capacityfor direct control of a vvide array of 
induslrial and aulomutive aclualcrs. Faoh rnodulB leaUires an inlegraled 
connector junction box with screvv-lerminal or cable options for flexible, 
luw-cusl signal wiring. All modules feature [he Nl CompaclRSQ Fxdreme 
Industrial Certifications and Ratings. 

System Compatibility 

You can use N! C Series modules interchangeably in CcmpactRIO and 
Nl CompactD.AQ chassis. Many of the advanced timing features 
described apply only Io CompaciRIO reconfigurable l/'O systems and not 
to Nl CampaclDAQ, 

Advanced Features 

When used in CcmpactRIO, C Series digital output modules connect 
directly lu rBcunfigyrablG I/O (RIO) FPGA hardware to create high- 
performance embedded systems. The reconfigurable FPG.A hardware 
within CcmpactRIO provides a variety of options for timing, triggering, 
synchronization, digital waveform generation, or digital communication. 
For instance, with CompaciRIO you can implement a circuit lo generate 


pulse-widtii mudulalion (PWM) outputs for controlling mclors, liealers, 
or tans, as well as perform pulse code modulation encoding (PCMF.) for 
wireless telemetry applications. 

Key Features 

• High-performance digital output switching for any CompactRiO embedded 
system, R Series e.xpansion chassis, or Nl CompactD.AQ system 
Screw'-terminal, strain relief, high-voltage, cable, solder-cup 
backsholl, and other connoctivity options 

• Channrd-to-eartli ground double-isolation barrier for safely 
and noise immunity 

• Nl CompactRiO Extreme Industrial Certifications and Ratings 

Visit iii.cofii/compactrw or for up-to-date 

information on module availability, example programs, application notes, 
and other developer tools. 

Typical certificatiuns - Actual specincations vary from product to 
product. Visit niconi/certificatbii for details. 
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Figure 86. Datasheet for C Series Digital Output and Relay Modules from ni.eom. 
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D. NI DEVELOPER SUITE WITH REAL-TIME AND FPGA SOFTWARE 

FOR COMPACTRIO 

From the CompactRIO ADVISOR on the webpage for this software paekage, 
http://sine.ni.eom/nips/cds/view/p/lang/en/nid/14294 “The NI Developer Suite Core 
Package with the Real-Time Deployment and FPGA Deployment options include all of 
the tools you need for building a custom embedded system - Lab VIEW FPGA for 
customizing hardware I/O in the user-configurable FPGA core and Lab VIEW Real-Time 
for building deterministic and reliable real-time applications. The NI Developer Suite 
Core Package with the Real-Time Deployment and EPGA Deployment options also 
include a comprehensive set of Lab VIEW add-on tools for algorithm development, real¬ 
time development analysis, Microsoft Office report generation, and database 
communication.” 

Of particular note in this software package is the EPGA Deployment Option, 
which includes the LabVIEW EPGA Module and NI RIO. These are essential 
programming components to the LabVIEW software library used with the ACS. NI RIO 
is needed to communicate with the CompactRIO Chassis(EPGAs) and the LabVIEW 
EPGA module for programming commands that can be auto-compiled into VHDL for the 
CompactRIO Chassis(PPGAs). 
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APPENDIX 1. LAB VIEW EOR THE DIGITAL PHASE SHIETER 


How the digital phase shifter is controlled by code written in Lab VIEW and run on the 
CRIO is explained here. 

The Lab VIEW Virtual Instruments (Vis) described here include: ManualControl, 
PhaseShifterControl, MasterControl, and SlaveControl. Together they provide three 
distinction functionalities in the development of the ACS’s software. These are the main 
Vis. All Vis, including Block Diagram and Eront Panel files for each, are included after 
the description of these main Vis. The Vis not described are test components only and 
not used with the ACS’s current functionalities, although they did play a part in its past 
development and possibly future development and are included here for that purpose. 
Primarily they were used to test certain functions that were later incorporated into the 
main Vis or demonstrate alternative ways of implementing the same functionality as was 
incorporated into the main Vis as part of testing method efficiency. These Vis include: 
datalO, fpgaRX_test, fpgaTX_test, NetworkRX_test, NetworkTX_test, and tDcalc. 
KNORR-CAL_project is an overall view of the Lab VIEW project with all Vis and 
hardware indicated. The VI PSC6-bit is merely an example of how to modify the files 
for a different number of control bits. When the ACS was first developed a 6-bit digital 
phase shifter rather than an 8-bit digital phase shifter was used. The basic principle of 
this modification is to globally find and replace the number representing 2^(current 
number of bits) with 2^(future number of bits). Eor example, replace 255 (for 8-bit) with 
63 (for 6-bit). 

Observe the Block Diagram for each VI and follow along in the diagram left to right as 
the description explains the graphical code. 

A. PHASESHIFTERCONTROL 

Eirst, the count is reset (initialized) to zero. Then, the time delay (tD) is read in to 
become the time between phase state changes. This is also referred to as ‘ticks’. 1 tick = 
1 bit change = 25 ns (the time of one 40MHz clock cycle). The count is checked so as to 
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not exceed the number of phase shifter bits (phase states). Here there are 8-bits or 
2^8=256 phase states. Since 0 is considered the first phase state, the count will increment 
through all 256 phase states and stop at 255. 

If the Doppler shift is going to be an up Doppler shift (increase in frequency) the 
count proceeds unaltered, counting up. If the Doppler shift is to be a down Doppler shift 
(decrease in frequency) the count proceeds inverted, counting down. The integer value of 
the count is then converted into a Boolean value. For example, the integer value of 157 = 
10011101. The Boolean value is then broken into each of its individual bits and routed to 
the corresponding Digital I/O (DIO) pin to apply or remove voltage from the appropriate 
phase shifter control pin. In this Vi’s current state, both DIO modules receive the same 
signal. 

See Chapter on Future Work. Once the programmable attenuator is received, the 
second group of DIO pins will be assigned to control this device with dependent 
statistics. For independent statistics, a second PhaseShifterControl VI will need to be 
created(duplicated and renamed AttenuatorControl). The only difference between the 
original and its duplicate is that the duplicate DIOs will be deleted one from each VI so 
that each VI contains only one unique set of DIO assignments. For dependent 
(correlated) statistics the DIOs will be controlled from within the same VI in a 
synchronized time delay (tD). 

Once the commands for the appropriate voltages have been sent, the count 
increments, and the process begins all over again. This process continues until one 
complete cycle through all the phase states is completed (one period of 2n). Then, the 
loop encounters an exit condition and the count is reset and the whole sequence begins 
again (the next period of 2n). 

B. SLAVECONTROL 

Paired with MasterControl provides automatic PhaseShifterControl to work with 

data read in from a file. The difference is that an array of integer values corresponding to 

the tD is first read and stored. Then, the individual time delay integer values are used to 

provide the appropriate phase shifter Doppler shift values as opposed to the manual entry 
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input by a user in the GUI of PhaseShifterControl. In addition, instead of a switeh being 
flipped by a user in the GUI to deeide up or down Doppler shift, the sign of the integer 
value is used. A negative sign is used to indieate a down shift and the absenee of a 
negative sign (i.e. positive) to indieate an up shift. Zero indieates no Doppler shift. The 
proeess is then exaetly the same as that implemented/programmed in PhaseShifterControl 
(See seetion above by the same name). 

*NOTE; the number of Doppler frequeneies eannot be automatieally deteeted from the 
data file stored on the PC. The number of frequeneies must be defined from within the 
VI. This allows the appropriate sized array to be ereated in the software even before it 
starts running. This is noted in the VI Front Panel (GUI) as well. 

C. MASTERCONTROL 

This VI reads in the value of the Doppler frequeneies from a tab delaminated file 
(sueh as how the array of variables would be exported/saved from the Weather Signal 
Simulator in MATLAB). The numbers(frequeneies) are read in as double floating point. 
numtDs is the number of Doppler frequeneies read in. The Doppler frequeneies (DF) are 
then used to ealeulate the Radial Veloeity (Vr), Period (T), and Bit Delay (t) aeeording to 
the following relations: 

Q 

Equation 1 Vr = DF x- in meters per seeond 

RF x2 

Equation 2 T = —x 1000 in milliseeonds 
DF 

T lO*” 

Equation 3 t = -in ‘tieks’ (of 25ns) 

n 25 

Where the RF Frequeney (RF) is defined to be 9990000000Hz or 9.99GHz, the Speed of 
Fight (e) is 299792458 m/s, and the number of Phase Steps (n) is 63 (for 6 bits). 

However, sinee only integers (INT) ean be passed to the FPGA, the doubles 
(DBF) are eonverted to integers and sent to the FPGA. The Period, DF, and Vr are then 
reealeulated based on these 16-bit integer values. This is to show the user what aetual 
values will be produeed by the ACS. This only affeets values stored in the data file 
whieh may have been truneated by the eonversion. This is a potential software limitation. 
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However, this will not affeet the ealibration process, only “real world” weather 
simulation, because any variation from the input values that have been truncated are now 
recomputed and displayed to the user. The radar should be tested/calibrated to agree with 
these values, not the values output by the Weather Signal Simulator. This is a limitation 
of the calibrator’s software when interacting with the FPGA. Although values may be 
calculated accurately by the Weather Signal Simulator to agree with “real world” weather 
patterns, the calibrator cannot always reproduce these exact values. This is a hardware 
limitation due to the FPGA and seen in the software. In practice, this limitation is similar 
to other hardware limitations already discussed. See Hardware Limitations Section. 

*NOTE: before executing the VI, the array size in SlaveControl must be defined to the 
number of Doppler frequencies intended to be read from the data file. This is not 
automatic. It cannot predetermine the number of entries. There is no “read until end of 
file.” See also SlaveControl NOTE. 

*NOTE; the path on the PC to the data file must be entered into the VI. The data file 
must be a single row of tab delaminated Doppler frequency entries stored in a text file 
(.txt) where a negative sign indicates a down Doppler shift and no sign indicates an up 
Doppler shift. Do not use + signs. 
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CODE 

Ul 


LabView Code (# indicates state change) 
User Interface 


datalO 


reads in data file containing Doppler Shift values input from 
Ulf Shroder's MATLAB weather simulator 


fpgaRX_test test file to send data read on laptop to fpga 

fpgaTX_test test file to receive data sent from laptop on fpga 

KNORR-CAL project overall software view. LabView Project. 

vi to manually control individual phase shifts. High or low 
any one or combination of bits 

vi to automatically control phase shifts. Runs on laptop, 
test file to send data read on laptop to CRIO 
test file to receive data sent from laptop on CRIO 
only difference from PhaseShifterControl is 6-bit (63) test 
PhaseShifterControl6-bit vs. 8-bit (255) 

fpga resident to control phaseshifter through a series of 
PhaseShifterControl phase shifts 

SlaveControl version of PhaseShifterControl to work with MasterControl. 

LabVIEW computer for the mathematical relations between 
time delay tD, radial velocity, and Doppler frequency. 
MATLAB Script Node is used in MasterControl instead of 

tDcalc this. 


ManualControl 
MasterControl 
NetworkRXtest 
NetworkTX test 


Table 10. Table Summary of Lab VIEW Vis. 
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Figure 87. datalO CODE 
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Figure 88. datalO UI 



Figure 89. fpgaRX test CODE 
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® fpgaRX test.vi Front Panel on FPGA Target 





Figure 90. fpgaRX test UI 



Figure 91. fpgaXX test CODE 



Figure 92. fpgaTX_test_UI 


136 



















































































































































































1 Project Explorer - KNORR-CAL.Ivproj 

i 


F8e Edit View Project Operate loots Window Help 

j| s 

X 



- 


B Prolect: KNORR-CAL.Ivpro1 

$ S CRBB 

I MasterControl.vi 
I I datalO.vi 
I fpgaTX_test.vi 
I h NetworkTXJ:est.vi 
I 1^ host.vi 
I Dependencies 
^ 'it BuikJ Specifications 
6 ^ CRIO-KNORR (169.254.0.8) 

$ H FPGATarget (RIOO::INSTR, cRIO-9101) 
i B Digital Line Input and Output 
MBS NI9401 

I h i* DIO0(NI9401) 
i i* DIO 1(NI 9401) 

I h i* DIO2(NI9401) 
h i* DIO3(NI9401) 

I h i* DIO4(NI9401) 
h i* DIOS (NI9401) 
i h i* DIO6(NI9401) 

' i* DI07 (NI 9401) 
i B S NI94012 

1 i* DIOO 2 (NI 9401 2/DIOO) 
h i* DIOl 2 (NI 9401 2/DIOl) 

I i* DI02 2 (NI 9401 2/DI02) 
h i* DI03 2 (NI 9401 2/DI03) 
i i* DI04 2 (NI 9401 2/DI04) 
h i* DI05 2 (NI 9401 2/DI05) 
i i* DI06 2 (NI 9401 2/DI06) 

- i* DI07 2 (NI 9401 2/DI07) 

1 40 MHz Onboard Cbck 

I h PhaseShifterControl.vi 
i 13^ PhaseShifterControl6-bit.vi 
i r SlaveControl.vi 
h isi ManualControl.vi 
i h fpgaRX_test.vi 
r ® NI9401(Stot2) 
h i* NI 9401 2 (Slot 3) 
h ^ Dependencies 
' it Build Specifications 
I ■ isi tDcaIc.vi 
h lii Net\vorkRX_test.vi 
i »!• target - single rate - variables.vi 
f ai support - acquire data.vi 
BCJ variables - network - RT.Ivlib 
i h'J^ Array 
!■ data 
stop 

i Dependencies 
' it Build Specifications 


Figure 93. KNORR-CAL_project 
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Figure 94. ManualControl CODE 
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Figure 95. ManualControl UI 



Figure 96. MasterControl CODE 
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Figure 97. MasterControl UI 



Figure 98. NetworkRX_test_CODE 
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Figure 99. NetworkRX_test_UI 



Figure 100. NetworkTX_test_CODE 
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Figure 101. PhaseShifterControl6-bit_CODEl 



Figure 102. PhaseShifterControl6-bit_CODE2 
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Figure 103. PhaseShifterControl6-bit_UI 



Figure 104. PhaseShifterControl CODEl 
143 
















































































































Figure 105. PhaseShifterControl_CODE2 



Figure 106. PhaseShifterControl UI 
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Figure 107. SlaveControl CODEO 



Figure 108. SlaveControl CODEl 
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Figure 109. SlaveControl_CODE2 



Figure 110. SlaveControl Ul 
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Figure 111. tDcalc CODE 
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Figure 112. tDcalc UI 
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